AD-A02O  685 


AH  INVESTIGATION  OF  THE  DISTRIBUTION  OF  RADIATION 
SCATTERED  BY  OPTICAL  SURFACES 

Roland  V.  Shack,  et  al 

Arizona  University 


Prepared  for: 

Space  and  Missile  Systems  Organization 


August  1975 


DISTRIBUTED  BY: 


NatiMal  Ttckiiol  lifonRatiM  Service 
U.  S.  DEPARTMENT  OF  COMMERCE 


AN  INVESTIGATION  OF  THE  DISTRIBUTION  OF  RADIATION 
SCATTERED  BY  OPTICAL  SURFACES 


Roland  V.  Shack 
Jftses  F.  Harvey 


Optical  Sciences  Center 
The  University  of  Arizona 
Tucson,  Arizona  SS721 


August  1975 

Final  Report  for  Period  October  1973  to  December  1974 


Approved  for  public  release;  distribution  unlimited. 


Prepared  fox  the 


Air  Force  Systems  Command 
Los  Angeles,  California  90009 


( s 1 r \ 

\r\v~'- 


Space  and  Missile  Systems  Organization  \\  \j 


\ Vat- 


Best  Avail  ao!a  Ob 


national  technical 
INFORMATION  SERVICE 

U S Dc-pur  Imer.t  of  Commerce 


This  report  has  been  reviewed  by  the  Information  Office  fcH)  ansi  is 
releasable  to  the  National  Technical  Information  Service  (NTIS).  At 
NHSf  it  will  be  available  to  the  General  public,  including  foreign 
nations . 


This  technical  report/has  reviewed  and  is  approved  for  publication. 



EDWARD1CT5IETZ 7Cz?trW&F 
Program  Manager 
Survivability  Directorate! 


FOR  THE  COMMANDER 

1/^'  \ j 

•ii’  (C  - 

Jam?  A . NbRB^O  LtCo^^OSAF 

director,  SurvivabilityjMplrectorate 


VA’^i 


J 


UNCLASSIFIED 

security  classification  of  this  pace  ibti**  o««»  Ent*»rf; 


REPORT  DOCUMENTATION  PAGE 


report  number 

TR  #75-236 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


i.  GOVT  ACCESSION  NO.!  »■  RECIPIENT'S  CATALOG  NUMBER 


4.  TITLE  !*nrt  Subtttlr) 

AN  INVESTIGATION  OF  THE  DISTRIBUTION  OF 
RADIATION  SCATTERED  BY  OPTICAL  SURFACES 


S.  TYPE  OF  REPORT  A PERIOO  COVEREO 

Final  Report  for  Period 
October  1*  73  to  Dec.  1974 


«.  PERFORMING  ORG.  REPORT  NUMBER 


7.  AUTHOR!*; 

Roland  V.  Shack 
James  E.  Harvey 


I.  CONTRACT  OR  GRANT  NUMBER!*; 


F04701-74-C-0083 


*.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Optical  Sciences  Center 
University  of  Arizona 
Tucson,  Arizona  85721 


10.  PROGRAM  ELEMENT  .’ROJECT,  T ASK 
AREA  4 WORK  UNIT  NUMBERS 


12.  REPORT  DATE 

August  1975 


II.  CONTROLLING  OFFICE  NAME  AND  ADORESS 

Space  and  Missile  Systems  Organization 
Air  Force  Systems  Command 

Los  Angeles,  California  90009  


4.  MONITORING  AGENCY  name  ft  AQORESSftf  dl/foront  from  Controlling  OIUco)  i 15.  SECURITY  CLASS,  {of  thio  ropotl) 


IS.  NUMBER  OF  RAGE* 


Unclassified 


ISa.  DECLASSIFICATION/ DOWNGRADING 
SCHEDULE 


l«.  DISTRIBUTION  STATEMENT  <ol  thit  AP*pw; 


Approved  for  public  release;  distribution  unlimited. 


It.  KEY  WORDS  (Cttiffnu*  on  to «w«*  oldo  H n «»«a**wy  tgonttly  by  block  mml  w, 

Surface  Scatter 

Diffuse  Reflectance 

Surface  Properties  ; 

Optical  Surfaces 
Diffraction  Theory 


SO  ABSTRACT  (Continuo  on  rtvww  •!«#•  II  nocoo<^t y mod  fdmfffjr  bf  block  mmbot) 

-The  direct  application  of  Fourier  transform  theory 
scatter  allows  the  well-known  techniques  of  linear 
applied  to  the  scattering  process.  For  surfaces  w 
this  approach  results  in  a linear,  shift- invariant 
does  not  change  shape  with  variations  in  the  angle 
tus  aj(d  technique  used  for  measuring  the  scattered 
variety  of  optical  surfaces  as  a function  of  angle 


to  the  problem  of  surface 
systems  theory  to  be 
ith  well-behaved  statistics, 
scattering  functicn  that 
of  incidence.  The  appara- 
light  distribution  from  a 
of  incidence  is  discussed, 


liwr 


FORM 

I JAN  71 


COITION  OF  I NOV  «*  It  OBSOLETE  , 
S/N  9|01-*M-4*01  f 


JNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  !*h«i  M*  BWM*M9 


UNCLASSIFIED 

-H-IJMITY  CLASSIFICATION  OP  THIS  PAOC.'Wi «n  P«C  tnlwrt) 


20.  Continued. 

— - and  the  results  of  experimental  measurements  are  presented  and  compared 
with  theoretical  predictions,. _ 


t 


I o 


UNCLASSIFIED 

-'leeuairv  classification  of  this  FAeer»n«n  v«>* «:«» 


FOREWORD 


This  is  the  final  report  for  the  investigation  of  radiation 
scattered  by  optical  surfaces  by  the  Optical  Sciences  Center,  University 
of  Arizona,  Tucson,  Arizona.  The  work  was  performed  under  A5r  Force 
Contract  FC4701-74-C-0083.  The  program  was  carried  out  under  the  aus- 
pices of  the  Space  and  Missile  Systems  Organization,  Air  Force  Systems 
Command,  Los  Angeles,  California. 

The  principal  investigator  for  the  Optical  Sciences  Center  was 
Roland  V.  Shack.  The  project  engineer  and  co-investigator  was  James  E. 


Harvey. 


TABLE  OF  CONTENTS 


Page 


I.  INTRODUCTION 1 

Previous  Studies  of  Surface  Scattering  Phenomena  1 

II.  GENERALIZED  SCALAR  DIFFRACTION  THEORY  6 

The  Diffracted  Wave  Field  as  a Superposition 

of  Plane  Waves  7 

Initial  Conditions  7 

The  Direct  Application  of  Fourier  Transform  Theory  . . 9 

The  Traits  fer  Function  of  Free  Space 10 

The  Diffracted  Wave  Field  as  a Superposition 

of  Spherical  Waves  14 

The  Point  Spread  Function  14 

Huygens'  Principle 15 

General  Rayleigh-Sommerfeld  Diffraction  Formula  ....  16 
Geometrical  Configurations  of  the  Observation  Space  ....  17 

The  Diffracted  Wave  Field  on  a Plane 17 

The  Diffracted  Wave  Field  on  a Hemi  sphere 19 

Aberrations  of  Diffracting  Systems  22 

Shift  Invariance  of  the  Diffracted  Wave  Field  28 

Summary  32 


III.  SURFACE  SCATTER  THEORY 35 

Surface  Scatter  Phenomena  as  a Diffraction  Process  ....  35 
Effect  of  the  Scattering  Surface  upon  the 

System  Pupil  Function  36 

Scattered  Wave  Field  on  a Hemisphere 38 

The  Effective  Spread  Function  of  a 

Scattering  System 39 

Scattering  from  Optical  Surfaces  .....  39 

Assumptions  Concerning  the  Statistical  Properties 

of  an  Optical  Surface  .....  41 

The  Effective  Transfer  Function  of  a 

Scattering  Surface  42 

The  Effective  Spread  Function  of  a 
Scattering  Surface  45 


ii 


♦ 


t 


TABLE  OF  CONTENTS — Continued 


9 


t 


I 


I 


Page 


Relationship  between  the  Spread  Function 

and  the  BRDF 47 

Shift  Invariance  in  Direction  Cosine  Space  48 

IV.  SURFACE  SCATTER  MEASUREMENTS  51 

Apparatus 51 

The  Light  Source 51 

The  Detector  Unit 54 

The  Scanning  Mechanism  60 

The  Incident  Beam 60 

Experimental  Procedures  60 

Sample  Preparation 64 

System  Alignment . 66 

Measurement  Technique  68 

Data  and  Results 71 

V.  SURFACE  STRUCTURE  MEASUREMENTS  78 

Electron  Microscopy  78 

FECO  Interferometry SJ 

VI.  SUMMARY 83 

APPENDIX  A.  FRESNEL  AND  FRAUNHOFER  DIFFRACTION 86 

APPENDIX  B.  CALCULATION  OF  ABERRATION  COEFFICIENTS  FOR 

DIFFRACTED  WAVE  FIELDS 89 

APPENDIX  C.  BIDIRECTIONAL  REFLECTANCE  DISTRIBUTION  FUNCTION  ...  94 

APPENDIX  D.  LITERATURE  REVIEW  97 


%J 


LIST  OF  ILLUSTRATIONS 


Figure  Page 

1.  Schematic  Representation  of  Scattering  from  a 

Rough  Surface 2 

2.  Geometry  of  Planes  P0  and  P 8 

3.  Unit  Circle  in  Direction  Cosine  Space  12 


4.  (a)  Original  Disturbance  Separated  into  its  Propagating 

and  Evanescent  Parts,  (b)  Direction  Cosine  Spectrum 
Separated  into  its  Propagating  and  Evanescent  Parts  ....  13 

5.  Geometrical  Relationship  between  Incident  Beam, 

Diffracting  Aperture,  and  Observation  Hemisphere  20 

6.  Diffraction  Pattern  of  a 10-Line  per  mn  Ronchi  Ruling 

Placed  in  an  f/6  Cone  of  Light  with  a 40-mm  Diameter  ...  27 

7.  Diffraction  Pattern  of  a 10-Line  per  mm  Ronchi  Ruling 

Placed  in  an  f/20  cone  of  Light  with  a 12-mra  Diameter  ...  29 


8.  Geometrical  Configuration  when  the  Incident  Beam 

Strikes  the  Diffraction  Aperture  at  an 

Arbitrary  Angle  31 

9.  Illustration  of  the  Position  of  the  Diffracted  Orders 

in  Real  Space  and  Direction  Cosine  Space  33 

10.  Graphical  Illustration  of  Surface  Height  Variations 

and  Associated  Statistical  Parameters  37 

11.  An  Ensemble  of  Two-Dimensional  Sample  Functions 

Representing  Surface  Height  Variations 
Constitutes  a Random  Process  W(x,y)  40 

12.  (a)  Effective  Transfer  Function  of  a Scattering  Surface. 

(b)  Spread  Function  Associated  with  a Scattering 

Surface 46 


IV 


LIST  OF  ILLUSTRATIONS- -Continued 


Figure  Page 

13.  Geometrical  Relationship  between  Incident  Beam, 

Scattering  Surfrce,  the  Measured  Scattered  Light 
Distribution,  and  the  Resulting  Spread  Function  49 

14.  Schematic  Diagram  of  Scatter  Measurement  Apparatus  52 

15.  Apparatus  for  Measuring  Scattered  Light  Distribution 

from  Optical  Surfaces  53 

16.  (aj  Previous  Configuration  of  Detector  Probe  Unit. 

(b)  New  Configuration  of  Detector  Probe  Unit 55 

17.  Detailed  Illustration  of  Hew  Fiber-Optic  Probe  57 

18.  Detector  Response  as  a Function  of  Field  Angle  for 

Detector  Configurations  Shown  in  Fig.  16  58 

19.  Photomultiplier  Tube  Linearity  Curve  59 

20.  Detailed  Illustration  of  Movable  Arm  Supporting  the 

Optical  System  Forming  Incident  Beam (>1 

21.  Geometrical  Configuration  of  Two  Principal  Planes 


in  which  the  Scattered  light  Field  will  be  Sampled  ....  63 


22.  (a)  Spread  Function  of  Incident  Beam,  (b)  Spread 

Function  of  Scattering  Surface,  (c)  Spread 

Function  of  Scattering  System  70 

23.  Comparison  of  Scatter  Profiles  for  Different  Incident 

Angles 72 

24.  Comparison  of  Scatter  Profiles  Taken  in  Different 

Directions  for  an  Incident  Beam  at  45  Degrees 74 

25.  Comparison  of  Scatter  Profiles  for  Samples-  with  a Wide 

Range  of  Surface  Roughnesses  for  a Normally  Incident 

Beam 76 

26.  Comparison  of  Scatter  Profiles  from  Special  Materials 

and  Unusual  Fabrication  Techniques  with  Those  of 
Conventional  Optical  Surfaces  . 11 


v 


LIST  OF  ILLUSTRATIONS — Continued 


Figure  Page 

27.  (a)  Typical  Electron-Micrograph  Stereo  Pair  of  a 

Ground  Glass  Surface,  (b)  Surface  Profile  of 

Ground  Glass  Sample  as  Determined  by  Stereo- 

Photogrammetric  Techniques  79 

28.  Surface  Profile  of  Ground  Glass  Sample  and  its 

Associated  Surface  Height  Distribution  and  j 

Autocorrelation  Function  81 

29.  Geometry  of  Incident  and  Reflected  Elementary  Beams 

used  to  Define  the  Bidirectional  Reflectance 

Distribution  Function  95 


vi 


?*r,&S£*^%$^%BSZ&!&P^3& 


LIST  OF  TABLES 


Table  Page 

1.  Tabulation  of  Expressions  for  the  Aberration 

Coefficients  for  Several  Different  Geometrical 

Configurations  of  (a)  Incident  Wavefront,  and 

(b)  Observation  Space  24 


2. 


Surface  Preparation  Techniques 


65 


CHAPTER  I 


INTRODUCTION 

Kith  the  recent  refinement  of  optical  devices  and  systems  there 
has  evolved  a need  to  examine  in  closer  detail  the  mechanism  and  proper- 
ties of  optical  scattering.  The  relationship  between  surface  micro- 
reughness  and  radiant  energy  scattering  has  an  important  bearing  on  the 
cost  and  performance  of  optical  systems.  The  surface  roughness  and 
resulting  light  scattering  properties  are  therefore  important  character- 
istics of  any  material  considered  for  use  as  an  optical  element.  More- 
over, if  the  scattering  mechanism  were  completely  understood,  surface 
preparation  techniques  could  possibly  be  changed  to  save  both  time  and 
expense. 


Previous  Studies  of  Surface  Scattering  Phenomena 
Tf  a propagating  wave  is  incident  upon  a plane  surface,  the 
reflected  wave  is  concentrated  in  the  specular  direction  as  determined 
by  the  well-known  laws  of  reflection  [Bom  and  Wolf  (1959)].  Another 
idealized  surface  is  the  perfectly  diffuse  reflector  which  scatters 
light  according  to  Lambert's  cosine  law  [Hudson  (1969)].  A more  physi- 
cally realistic  situation  is  shown  in  Fig.  1,  which  illustrates  the 
optical  scattering  that  occurs  when  light  is  reflected  from  a rough  sur- 
face. If  the  surface  is  not  too  rough  the  reflected  light  is  seen  to 
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consist  of  a specular  component  and  a diffuse  component  which  is  scat- 
tered over  a wide  range  of  angles  centered  upon  the  specular  beam. 

One  of  the  earliest  investigators  of  scattering  from  n rough 
surface  was  Lord  Rayleigh.  In  1896  [Rayleigh  (1945) J he  was  investiga- 
ting the  reflection  of  acoustic  waves,  and  later  [Rayleigh  (1901) J he 
noted  the  effects  of  poorly  polished  surfaces  on  optical  performance. 

His  work  examined  the  effects  of  surface  roughness,  wavelength,  and 
angle  of  incidence  on  the  reflected  beam.  Chenmoganadam  (1919)  derived 
a theory  of  scattered  light  based  on  the  phase  shift  of  the  reflected 
beam  due  to  the  rough  surface. 

However,  it  was  not  until  the  problem  of  background  clutter  in 
radar  applications  became  apparent  that  a determined  effort  was  made  to 
solve  the  scattering  problem  for  random  surfaces.  For  example,  scatter- 
ing from  the  sea  motivated  the  work  of  Davies  (1954)  as  well  as  others 
[Blake  (1950),  Barrick  (1970),  Bass  (1968),  Beard  (1961),  and  Fuks 
(1966)].  Considerable  work  has  also  been  done  in  attempts  to  explain 
radar  reflection  from  the  mocn  [Daniels  (1961),  Evans  and  Pcttengill 
(1963),  Fung  (1964,  1967),  liagfors  (1964)]. 

Random  rough  surfaces  can  be  classified  into  two  main  groups. 
Rough  surfaces  made  up  of  a random  array  of  objects  or  shapes  with  known 
scattering  characteristics  were  investigated  by  Ament  (1960),  Twersky 
(1957),  Spetner  (1958),  and  Peake  (1959).  The  other  approach  treats  the 
rough  surface  as  a stochastic  process.  This  approach  was  taken  by 
Isakovich  (1952),  Ament  (1953),  Eckart  (1953),  Feinstein  (1954),  Davies 
(1954),  and  Beckmann  (19S7). 
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Since  optical  surfaces  clearly  fall  into  the  second  classifica- 
tion of  randoa  surfaces,  Bennett  and  Porteus  (1961)  expanded  and 
experimental ly  investigated  the  scattering  theory  of  Davies  (1954) . 

From  this  and  subsequent  work  [Bennett  (1963)  and  Porteus  (1963)]  the 
reflectance  properties  of  samples  with  a measured  surface  roughness  were 
directly  compared  to  theory  with  good  results.  Interest  in  these  mea- 
surements led  to  investigations  at  the  Optical  Sciences  Center  by 
McKenney  (1972),  Mott  (1971),  Orme  (1972)  and  DeBell  (1974). 

In  the  above  studies,  satisfactory  theories  have  been  developed 
only  for  the  two  limiting  cases  of  very  rough  surfaces  (o/X  » 1)  and 
very  smooth  surfaces  (o/X  « 1).  Until  recently  most  of  the  efforts 
have  been  confined  to  very  smooth  optical  surfaces  which  satisfy  the 
second  condition.  However,  recent  applications  require  a more  general 
theoretical  treatment. 

In  this  study  surface  scatter  phenomena  is  treated  as  a diffrac- 
tion process  where  the  scattering  surface  introduces  an  isotropic  random 
wavefront  disturbance.  The  Fourier  techniques  of  linear  systems  theory, 
which  have  been  applied  to  the  problem  of  image  formation  through  a 
turbulent  atmosphere  by  Shack  (1967),  are  then  utilized  to  develop  a 
surface  scattering  theory.  Experimental  measurements  are  also  made  in 
an  attempt  to  verify  this  theory. 

Early  work  in  the  development  of  the  above  theory  was  carried 
out  under  Contract  F04701-72-C-0181  and  is  documented  along  with 
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expei iraental  results  in  "Surface  Scatter  Study,"  Technical  Report 

SAMSO  TR  74-8. , dated  February  1974.  An  update  of  the  literature  review 

accomplished  under  that  contract  is  included  in  this  report. 
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CHAPTER  II 


GENERALIZED  SCALAR  DIFFRACTION  THEORY 


The  phenomenon  of  diffraction  involves  a wave  field  incident 
upon  one  or  more  obstacles  or  apertures  of  absorbing  or  conducting  sur- 
faces. The  calculation  of  the  wave  field  emerging  from  such  a dif- 
fracting system  is  the  goal  of  all  diffraction  theories. 

It  should  be  emphasized  that  both  the  Kirchhoff  and  Payleigh- 
Sommerfeld  theories,  as  well  as  the  present  discussion  in  this  paper, 
treats  light  as  a scalar  phenomenon.  (For  a detailed  treatment  of  the 
historical  development  of  diffraction  theory,  see  Goodman  (1968), 
pp.  30-56.)  Such  an  approach  entirely  neglects  the  fact  that  the  vari- 
ous components  of  the  electric  ?Jid  magnetic  field  vectors  are  coupled 
through  Maxwell's  equations  and  cannot  be  treated  independently. 

Microwave  experiments  have  shown  that  scalar  theory  yields  very 
accurate  results  provided  that:  (i)  the  diffracting  aperture  is  larg* 

compared  to  a wavelength,  and  (ii)  the  diffracted  wave  field  is  ol served 
far  from  the  aperture.  It  is  significant  that  although  die  present 
treatment  is  limited  by  being  a scalar  theory,  the  above  approximations 
need  not  be  imposed  during  the  mathematical  formulation  as  they  are  in 
the  Kirchhoff  theory,  or  in  its  subsequent  use  as  they  almost  always  are- 
with  the  Rayleigh-Smamerfeld  theory.  The  following  development 
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threfore  provides  much  more  insight  and  accuracy  than  the  conventional 
treatments . 


The  Diffracted  Wave  Field 
as  a Superposition  of  Plane  Waves 

The  fundamental  diffraction  problem  basically  consists  of  two 
parts:  (i)  determining  the  effect  of  introducing  the  diffracting  screen 

upon  the  field  immediately  behind  the  diffracting  screen,  and  (ii)  how 
does  it  affect  the  field  downstream  from  the  diffracting  screen  (i.e., 
what  is  the  field  immediately  behind  the  diffracting  screen  and  how  does 
it  propagate}. 

Consider  first  the  propagation  problem  and  let  the  complex 
amplitude  distribution  of  the  optical  disturbance  in  plane  PQ  be  repre- 
sented by  the  scalar  function  U0(x,y; 0).  This  scalar  disturbance  in  P0 
will  be  considered  the  only  radiation  contributing  to  the  field 
U{x,y;z)  in  plane  P (see  Fig.  2).  z has  a parametric  relationship  since 
it  is  a function  of  the  observation  plane.  Note  that  a scaled  coordi- 
nate system  is  utilized  in  which  x = x/A,  y - y/\,  z = a/ A. 

Initial  Conditions 

It  will  be  assumed  that  the  complex  amplitude  of  any  monochroma- 
tic optical  disturbance  propagating  through  free  space  must  obey  the 
time- independent  wave  equation  (Helmholtz  equation).  We  will  also 
assume  that  the  Fourier  transform  of  the  scalar  field  U0(x,y; 0)  exists. 
This  is  not  a severe  restriction,  however,  as  Bracewell  (1965)  points 
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uo(*.-y;0) 

Fig.  2 


U(x,y;z) 


. Geometry  of  Planes  Pq  and  P. 
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out  that  physical  possibility  is  a valid  sufficient  condition  for  the 
existence  of  a Fourier  transform. 


The  Direct  Application  of  Fourier  Transform  Theory 

We  can  thus  define  the  following  Fourier  transform  relationships 
that  exist  for  planes  PQ  and  P. 


do(a,8;0)  = ||  Vo{x,y\0)e~i2*^*^dxdy 

—oo 

oo 

Ua(x,y; 0)  = ||  Ao(a,B;0)e£2w(fl*+6£WB 

— 00 

Od 

A(a,8;S)  = ||  U(x,yiz)c~*2*^+^dxdy 


U(x,y;z)  = ||  4(a,8;a)e^2*^+6^'£fa<ig. 


(1) 


(2) 


t3) 


(4) 


Equations  (2)  and  (4)  indicate  that  the  monochromatic  scalar  wave 
field  in  planes  P0  and  P can  be  decomposed  into  plane  wave  components 
whose  amplitudes  are  a function  of  the  direction  cosines  of  the  propaga- 
tion vector.  The  functions  Ao(a,6;0)  and  A(a,8;z)  will  be  referred  to 
as  the  direction  cosine  spectrum  of  plane  waves  contributing  to  the  dis- 
turbance U0(x,y; 0)  and  U(x,y’,z ) respectively.  The  direction  cosine 
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spectna  of  plane  waves  is  used  here  in  lieu  of  the  angular  .spectrum  of 
plane  waves  discussed  by  Ratcliff  (1956)  and  others.  This  is  consistent 
with  a more  general  treatment  which  is  not  restricted  to  small  angles. 

A A 

In  the  scaled  coordinate  system  V2  * X2V2,  and  k2  = \2k2  = (2ir)2. 
Hence  the  Helmholtz  equation  becomes 

[72  + (2  «)2]U{x,y',z)  = 0.  (5) 

Now  applying  Ea.  (4)  and  requiring  the  individual  plane  wave  components 
to  satisfy  the  Helmholtz  equation*  we  find 

yKct.Bja)  = Ao:  a,&;0)ei2vyz  (6) 

where 

Y = A - a sr-r-gr. 

The  Trans "er  Function  of  Free  Space 

Since  Eq.  (6)  relates  the  Fourier  transforms  of  the  scalar 
fields  in  planes  Pa  and  P it  can  be  rewritten  in  terms  of  a transfer 
function  for  free  space,  a,8;S) 


H(a,6;z) 


i4(a,g;a)  Unyz 

Ao( a,g;Q)  = e 


(7) 


We  have  thus  far  applied  no  restrictions  on  y and  two  regions  of 
interest  are  apparent:  that  for  real  values  of  y i-'nd  that  for  imaginary 
values. 
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for  (a2  + 32)  s 1 

1 


for  (a2  ♦ $2)  > 1 


Y is  real 

Y is  imaginary. 

(8) 


Consider  now  a unit  circle  in  the  a- 3 plane  of  direction  cosine 
space  as  shown  in  Fig.  3.  Inside  this  unit  circle  y is  real  and  the 
corresponding  part  of  the  disturbance  will  propagate  and  contribute  to 
the  wave  field  in  plane  P.  However,  those  components  of  the  direction 
cosir.e  spectrum  which  lie  outside  the  unit  circle  have  imaginary  values 
of  y and  represent  that  part  of  the  disturbance  which  experiences  a 
rapid  exponential  decay.  This  is  the  part  of  the  disturbance  which  is 
commonly  referred  to  as  the  evanescent  wave  [Goodman  (1968)]. 

Let  U0[x,yi 0)  be  the  product  of  the  complex  amplitude  transmit- 
tance of  a diffracting  screen  and  the  complex  amplitude  distribution 
incident  upon  it.  Figure  4a  illustrates  this  quantity  broken  down  into 
the  part  which  propagates  and  the  part  which  makes  up  the  evanescent 
wave  for  the  case  of  a unit  amplitude  j..’!we  wave  normally  incident  upon 
a circular  aperture.  The  direction  cosine  spectrum  of  plane  waves 
associated  with  these  respective  optical  disturbances  are  shown  in 
Fig.  4b. 

Note  that  the  sharp  corners  on  the  original  disturbance  in 
Fig.  4a  correspond  to  Kirchhoff's  unwarranted  boundary  conditions.  It 
is  the  propagating  part  only  that  accurately  represents  the  disturbance 
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Fig  3.  Unit  Circle  in  Direction  Cosine  Space 

The  plane  wave  components  inside  this 
circle  will  propagate,  and  the  plane 
wave  components  outside  this  circle 
contribute  to  the  evanescent  wave. 
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immediately  behind  the  diffracting  aperture  which  will  contrioute  to  the 
disturbance  downstream. 

It  is  now  clear  that  the  complex  amplitude  distribution  in  plane 
P can  be  determined  by  Fourier  transforming  the  original  disturbance 
Uo(x,y;0),  then  multiplying  the  resulting  direction  cosine  spectrum  of 
plane  waves  i4o(a,g;0)  by  the  transfer  function  of  free  space  given  in 
Eq.  (7),  and  finally  by  applying  the  inverse  Fourier  transform  integral 
of  Eq.  (4).  However,  the  limits  of  integration  on  Eq.  (4)  must  be 
changed  such  that  the  integration  is  performed  only  over  the  unit  circle 
instead  of  over  the  entire  a-B  plane. 

The  above  analysis,  in  which  an  optical  disturbance  is  repre- 
sented as  a superposition  of  plane  waves,  corresponds  to  the  transfer 
function  apprtach  in  image  formation  and  yields  considerable  insight 
into  the  behavior  of  these  plane  wave  components  during  the  phenomenon 
of  diffraction. 


The  Diffracted  Wave  Field 
as  a Superposition  of  Spherical  Waves 
The  convolution  theorem  [Bracewell  (1965)]  of  Fourier  transform 
theory  requires  that  a convolution  operation  exists  in  the  domain  of 
real  space  that  is  equivalent  to  Eq.  (6). 

The  Point  Spread  Function 

We  thus  have  the  alternative  method  of  expressing  the  complex 
amplitude  distribution  in  the  observation  plane  by  the  convolution  of 
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the  original  disturbance  with  a point  spread  function.  The  point  spread 
function  is  obtained  by  taking  the  inverse  Fourier  transform  of  the 
transfer  function  found  in  Eq.  (7). 

Starting  with  the  well-known  Weyl  expansion  formula 
[Weyl  (1919)],  Lalor  (1968)  obtained  a result  which,  with  straightfor- 
ward modification,  yields 


-2ir  j|  ei2,(ait^«, 


(9) 


where 


?2  a x2-  + y2  + z2. 


Hence  the  appropriate  point  spread  function  is  giver,  by 


Hx,y-,z)  = JT-I  jei2irYSj  = (J-.Ai 

\2tti*  / r 


i2*r 


(10) 


Huygens ' Princip le 

Recall  now  the  assertion  by  Christiaan  Huygens  [Thompson  (1912)] 
in  1678  that  each  element  of  a wavefront  may  be  regarded  as  the  center 
of  a secondary  disturbance  which  gives  rise  to  spherical  wavelets;  and 
moreover  that  the  position  of  the  wavefront  at  any  later  time  is  the 
envelope  of  all  such  wavelets.  These  intuitive  convictions,  sometimes 
called  Huygens'  wavefront  construction,  are  an  excellent  description  of 
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a convolution  operation  in  which  the  initial  disturbance  is  convolved 
with  a Huygens*  wavelet.  It  is  therefore  quite  appropriate  to  think  of 
the  point  spread  function  of  a diffraction  system  as  the  intersection  of 
a Huygens'  wavelet  with  the  observation  plane. 

Equation  (10)  is  therefore  an  exact  mathematical  expression  for  a 
Huygens'  wavelet  which  is  valid  right  down  to  the  initial  distuibance 
itself.  However,  for  r » 1, 

£2irr  £2ir(r  - 1/4) 

h(x,y\z)  = -t(z/r)  $ = (z/r)  , (II) 

r 

it  reduces  to  the  familiar  expression  for  a spherical  wave  with  cosine 
obliquity  factor  and  a u/2  phase  delay. 

General  Rayleigh-Sommerfeld  Diffraction  Formula 

If  we  write  down  th**  convolution  integral  for  the  disturbance  in 
the  observation  plane,  using  the  expression  in  Eq.  (10)  for  h(x,y;z ),  we 
obtain  the  general  Rsvleigh-Sommerfeld  diffraction  formula 

U(x,y;z')  = ||  IJ0{x',y'  ;0)  z/l  ^-7-^  dx'dy'  (12) 

— 00 

where 

i2  - (x-x1)2  + (y-y ')2  + z2.  (13) 

This  is  an  exact  expression  for  tne  diffracted  wave  field  which  is  valid 
throughout  the  entire  space  in  which  the  diffraction  occurs — right  down 
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to  the  aperture.  No  approximations  have  been  made.  Furthermore,  the 
above  equation  expresses  the  disturbance  on  the  observation  plane  as  a 
superposition  of  spherical  waves  which  corresponds  to  the  spread  func- 
tion approach  in  image  formation. 

Geometrical  Configurations  of  the  Observation  Space 
In  order  to  insure  a space  invariant  point  spread  function  our 
equations  have  been  restricted  to  mapping  an  optical  disturbance  from  an 
input  plane  to  an  output  plane,  where  £ has  a parametric  relationship 
since  it  is  a function  of  the  output  plane.  However,  the  summation  of 
these  Huygens'  wavelets  is  valid  over  anjl  surface.  The  above  treatment 
thus  gives  us  a far  more  powerful  conceptual  tool  than  provided  by  the 
equations  themselves. 

We  will  therefore  investigate  the  properties  of  the  diffracted 
wave  field  on  two  particular  geometrical  configurations  of  the  observa- 
tion surface. 

The  Diffracted  Wave  Field  on  a Plane 

Equation  (12)  reduces  to  the  more  familiar  but  less  general  form 
of  the  Rayleigh-Sommerfeld  diffraction  formula  [Goodman  (1968)]  when 
z » 1.  If  we  then  make  the  following  algebraic  substitutions 

I - 2(1  ♦ 6),  6 = , (14) 

z 

Eq.  (12)  can  be  rewritten  as 
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(15) 


# 00 
z2ir  2 
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Note  that  we  have  imposed  no  restriction  upon  the  size  of  the  aperture 
or  the  size  of  the  observation  space.  The  only  limitation  on  the  above 
equation  is  that  the  observation  plane  must  be  many  wavelengths  from  the 
aperture. 

The  above  diffraction  formula  is  a rather  unwieldy  integral  to 
solve  explicitly  for  most  problems  of  practical  interest.  Diffraction 
pattern  calculations  aro  greatly  simplified  by  using  Che  Fresnel  or 
Fraunhofer  diffraction  formula  (see  Appendix  A);  however,  severe 
restrictions  are  then  imposed  upon  the  size  of  the  aperture  and  the 
region  over  which  the  calculations  are  valid  in  the  observation  plane. 

In  order  that  we  do  not  impose  these  restrictions,  let  us  make  a 
binomial  expansion  of  the  quantity  Z in  the  exponent  of  Eq.  (15).  We 
can  then  rewrite  Eq.  (15)  as  the  following  Fourier  transform  integral 


• a 00 

z2  tiz 


Xltz  rr 

V{x,y\z ) = — — \\ty0{x' ,y' \x,y) 

iz  i i 


- % (*'*  + y'y) 


dx'dy' , 
(16) 


where  the  coop  lex  quantity 


3/0(*',$';x,P)  = To{x',y",0)  JI±jTei2vW, 


(17) 


can  be  regarded  as  a generalized  pupil  function.  A1 1 of  the  terms  from 
the  binomial  expansion  for  the  quantity  Z,  except  for  the  term  which  was 
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extracted  for  use  as  the  Fourier  kernel,  are  lumped  together  in  the 
quantity  ft  along  with  any  phase  variations  in  the  incident  wavefront. 

Equation  (16)  clearly  reduces  to  the  conventional  Fresnel  dif- 
fraction formula  when  a plane  wave  is  incident  upon  the  aperture  and 
when  z is  sufficiently  large  such  that  l is  adequately  approximated  by 
retaining  only  the  first  two  terms  of  the  binomial  expansion. 

The  Diffracted  Nave  Field  on  a Hemisphere 

Let  us  now  examine  the  diffracted  wave  field  on  a hemisphere 
centered  upon  the  diffracting  aperture  as  illustrated  in  Fig.  5.  The 
position  of  an  arbitrary  observation  point  will  be  specified  by  the 
direction  cosines  a and  B of  its  position  vector,  and  the  radius  r of 
the  hemisphere  upon  which  it  resides.  Note  that 

a ■ x/>,  6 = p/?,  and  y = a/F,  (18) 

where 

r2  = x2  + p2  ♦ a2.  (19) 

With  the  following  algebraic  substitutions 

l = r(l+c);  e = (£-r)/r.  (20) 
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Equation  (12)  can  be  rewritten  as 


• ^ <*» 

%2nr 


t/(o.g;r)  ° y 


2irr(l+e) 


- t 


(21) 


We  now  have  an  exact  expression  for  the  diffracted  wave  field  on  an 
observation  hemisphere  which  is  valid  throughout  the  entire  half-space 
behind  the  plane  of  the  diffracting  aperture. 

If  we  now  require  that  r » 1 and  make  the  appropriate  binomial 
expansions  for  the  quantity  (£-r) , we  again  obtain  a Fourier  transform 
integral 


00 

i2vr 


''(o.B;*)  = y 


ir 


\\%&',y';ai,Ve-i2*iax'  + ^^dx'ckj', 

(22) 


where  the  generalized  pupil  function  is  r,  ven  by 


«»0)  = T0&',y",G) 


(23) 


Once  again,  all  of  the  terms  from  the  binomial  expansion  for  the  quan- 
tity (£-r) , except  for  the  term  which  was  extracted  for  use  as  the 
Fourier  kernel,  are  lumped  together  in  the  quantity  W alorg  with  any 
phase  variations  in  the  incident  wavefront. 
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Aberrations  of  Diffracting  Systems 

The  quantity  W in  Eqs.  (17)  and  (23)  represent  phase  variations 
in  the  diffracted  wavefront  emerging  from  the  aperture.  Therefore,  fr 
can  be  interpreted  as  a conventional  wavefront  aberration  function 
[Hopkins  (1950)],  which  is  conveniently  expressed  as  a general  power 
series  expansion  of  the  pupil  coordinates  and  the  appropriate  field 
parameters . 

For  the  case  of  a rotationally-s>mmetric  diffracting  aperture  we 
can,  without  loss  of  generality,  choose  the  observation  point  to  lie  on 
the  y- axis  (x  = 0).  The  wavefront  aberration  function  can  then  be 
written  as 

* ^2Q0P2  + fro20«2  + frlllP®  cos$ 

+ frtOOP4  + ^040®4+  fc'niP^COS^  + &222P2<*2  COS2(J> 

+ ^220P2“  + ^31lP  “ COS 

+ higher-order  terms,  (24) 

where  p is  a normalized  field  position  of  the  observation  point  and  a 
is  a normalized  pupil  height. 

By  equating  coefficients  of  the  corresponding  terms  in  the 
appropriate  binomial  expansions  and  the  above  wavefront  aberration  func- 
tion, we  obtain  expressions  for  the  aberration  coefficients  in  terms  of 
the  aperture  diameter,  the  observation  distance,  and  the  appropriate 
field  parameter. 
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fhe  expressions  for  these  aberration  coefficients  are  derived  in 

Appendix  B and  tabulated  in  Table  1 to  enable  easy  comparison  of  several 

different  geometrical  configurations  of  the  incident  wavefront  and  the 
O 

observation  space. 

Consider  first  a plane  wave  illuminating  the  diffracting  aper- 
ture and  a plane  observation  space.  We  s^e  from  the  first  column  of 

O 

Table  1 that  all  aberrations  are  present  except  for  lateral  magnifica- 
tion error  (J?m)  which  is  absent  for  all  geometrical  configurations 

because  this  term  of  the  binomial  expansion  for  t is  extracted  for  use 

G 

as  the  kernel  of  the  Fourier  transform  integral.  It  is  clear  that  very 
large  observation  distances  are  required  to  reduce  defocus  (Vo 20)  to  a 
negligible  value.  Also,  distortion  (^311)  imposes  severe  restraints 
upon  the  field  angle  over  which  the  diffracted  wave  field  is  accurately 
described  by  the  Fourier  transform  of  the  aperture  function.  These 
restrictions  are  the  same  as  those  usually  imposed  during  the  develop- 
ment of  the  Kirchhoff  theory  and  in  most  applications  of  the  Rayleigh- 
Sommerfeld  diffraction  theory. 

The  effect  of  illuminating  the  aperture  with  a spherical  wave 
converging  to  the  observation  plane  is  to  eliminate  defocus  (Vo c 0 ) 
all  orders  of  spherical  aberration  (Voi*o) - This  removes  the  requirement 
for  an  extremely  large  observation  distance,  but  the  Fourier  transform 
of  the  aperture  function  is  still  valid  only  over  a small  region  about 
the  optical  axis  in  the  observation  plane. 
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Choosing  the  observation  space  to  be  a hemisphere  centered  upon 
the  diffracting  aperture  eliminates  field  curvature  (*220) » distortion 

A 

(*3ll)»  and  all  orders  of  piston  error  (*200  and  *400). 

Hence  for  the  case  of  a spherical  incident  wave  converging  to 
the  intersection  of  the  optical  axis  and  an  observation  hemisphere,  only 
coma  (&131)  and  astigmatism  (*222)  are  present.  And  the  value;  >f  the 
aberration  coefficients  can  be  calculated  from  the  relationships  pro- 
vided in  the  last  coluan  of  Table  1. 

Titus,  for  a system  with  an  aperture  diameter  of  1 mm  and  an 
observation  hemisphere  with  a radius  of  1 m,  we  have  for  X = 0.5  pm  and 

®max  “ 1 

V131  * 1.25  x 10“4 

£222  -2.50  x 10" l.  (25) 

Hence  there  is  only  X/4  of  astigmatism  at  the  edge  of  the  field  (i.e., 
90°  field  angle). 

A A 

Similarly  for  an  f/10  system  (d/r  = 0.1)  with  an  aperture  diam- 
eter of  5 cm  and  a maximum  field  angle  of  0.5°  (Bmax  “ 0.00872),  we  have 
for  X = 0.5  pm. 


*131  = 

0.545 

*222  * 

-0.095. 

(26) 

Hence  coma  dominates  at  the  edge  of  the  field  with  a value  of  approxi- 
mately X/2. 
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Finally,  for  an  f/10  system  with  an  aperture  diameter  of  1 cm 
and  a field  size  equal  to  the  size  of  the  aperture  (i.e.,  a pair  of  f/iO 
relay  lenses  1 cm  in  diameter),  we  have  for  X - 0.5  urn, 

i/l3l  = 0.625 

V222  = -0.625.  (27) 

We  find  slightly  more  than  A/2  of  coma  and  astigmatism  at  the  edge  of 
the  field. 

It  should  "be  pointed  out  that  in  each  of  the  above  cases  the 

radius  of  the  Airy  disc  in  direction  cosine  space  is  approximately  equal 

to  0.001  B . Hence  the  off-axis  aberrations  are  of  little  more  than 

max 

academic  interest  unless  there  is  some  structure  in  the  aperture  with 
high  spatial  frequency  content  which  will  diffract  light  at  large  angles 
from  the  direction  of  the  incident  beam. 

However,  the  above  analysis  of  the  aberrations  associated  with 
the  diffraction  process  can  be  readily  applied  to  holographic  systems  or 
systems  containing  diffraction  gratings. 

For  example,  an  f/6  system  with  a 10-line-per-mm  Ronchi  ruling 
placed  in  a 40-mrn  diameter  aperture  produces  the  diffraction  pattern 
shown  in  Fig.  6.  The  diffracted  order  at  three  different  field  posi- 
tions was  photographed  through  a microscope  with  the  following  results: 
at  B = 0 no  aberrations  were  apparent;  at  B = 0.04  coma  was  predominate 
with  a value  of  approximately  5X;  and  at  B - 0.10  coma  and  astigmatism 
both  have  values  of  approximately  ISA. 
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Fig.  6.  Diffraction  Pattern  of  a 10-Line  per  mm  Ronchi  Ruling  Placed  in 
an  f/6  Cone  of  Light  with  a 40-mm  Diameter. 

Magnified  images  of  diffracted  orders  at  various  field  positions 
indicate  that  coma  is  predominant  for  small  field  angles  with 
astigmatism  also  becoming  significant  at  larger  field  angles. 
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By  stopping  the  aperture  down  to  a diameter  of  12  mm  (f/20 
system)  and  observing  the  diffracted  orders  at  larger  field  angles, 
astigmatism  becomes  the  predominate  aberration  as  shown  in  Fig.  7.  The 
sagittal  focus  lies  on  the  observation  hemisphere  of  radius  * and  the 
medial  and  tangential  surfaces  have  smaller  radii  as  indicated.  The 
diffracted  order  at  M .020  exhibits  ?bout  6A  of  astigmatism.  This 
order  was  observed  through  focus  with  the  microscope  and  the  magnified 
images  are  displayed. 

In  both  of  the  above  examples  the  observed  aberrations  are  in 
good  agreement  with  those  predicted  from  the  coefficients  presented  in 
Table  1. 


Shift  Invariance  of  the  Diffracted  Wave  Field 
We  have  shown  that  any  departures  of  the  actual  diffracted  wave 
field  from  that  predicted  by  the  simple  Fourier  transform  of  the  aper- 
ture function  take  the  form  of  conventional  wavefront  aberrations. 

If  we  neglect  these  aberrations,  Eq.  (22)  reduces  to 


V{ o.Bjr)  = y ff  To(x',y';0)  e"^ir (ox'-S 

iv  ** 


(28) 


But  this  is  merely  the  Fourier  transform  of  the  aperture  function  multi- 
plied by  a spherical  Huygens’  wavelet. 


i2vr 

£/(a,B;r)  = y — — - — 3f{To(x,y',0)}. 

ir 


(23) 
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Fip.  7.  Diffraction  Pattern  of  a 10-Line  per  mm  Ronchi  Ruling  Placed 
in  an  i/20  Core  of  Light  with  a 12-mm  Diameter. 

Magnified  images  of  a diffracted  order  at  different  focal 
positions  indicate  that  astigmatism  is  predominant.  The 
relationship  of  Che  sagittal,  medial,  and  tangential  sur- 
face to  the  observation  hemisphere  is  also  shown. 


This  relationship  is  valid  provided  that  the  observation  space  is  a 
hemisphere  centered  on  the  diffracting  aperture  and  if  the  incident 
radiation  is  a unit  amplitude  spherical  wave  whose  center  of  curvature 
lies  on  the  intersection  of  the  observation  hemisphere  with  the  2-axis. 
Furthermore,  if  r is  large  compared  to  the  size  of  the  diffracting  aper- 
ture, the  Fourier  transform  relationship  is  accurate,  not  merely  over  a 
small  region  about  the  2-axis,  but  instead  over  the  entire  hemisphere. 

Now  consider  the  situation  where  the  incident  radiation  strikes 
the  diffracting  aperture  at  an  angle  6#  as  illustrated  in  Fig.  8.  This 
is  equivalent  to  introducing  a linear  phase  variation  across  the  aper- 
ture. By  applying  the  shift  theorem  [Bracewell  (1965)]  of  Fourier 
transform  theory  to  Eq.  (29)  we  find  that  the  complex  amplitude  distri- 
bution on  the  hemisphere  is  a function  of  (8  - Bo) » 

£2irr 

U(a*B-Boir)  = Y - — ~^{T(x}y;0)  exp[i2ir0o£] },  (30) 

ir 

where  3 is  the  direction  cosine  of  the  position  vector  of  the  observa- 
tion point,  and  80  is  the  direction  cosine  of  the  position  vector  of  the 
undiffracted  beam.  Note  that  these  direction  cosines  are  obtained  by 
merely  projecting  the  respective  points  on  the  hemisphere  back  on  to  the 
plane  of  the  aperture  and  normalizing  to  a unit  radius.  The  complex 
amplitude  distribution  at  an  arbitrary  point  on  the  hemisphere  can  now 
be  said  to  be  a function  of  the  distance  of  the  observation  point  from 
the  undiffracted  beam  in  direction  cosine  space . 
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Fig.  3.  Geometrical  Configuration  when  the  Indicent 
Beam  Strikes  the  Diffraction  Aperture  at  an 
Arbitrary  Angle. 
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As  a specific  example,  suppose  we  have  incident  light  striking  a 
diffraction  grating  at  an  angle  Qo-  The  diffracted  orders  will  strike 
the  observation  hemisphere  in  a cross  section  which  is  not  a great  cir- 
cle but  instead  a latitude  slice  as  shown  in  Fig.  9.  Thus  for  large 
angles  of  incidence  the  various  orders  appear  to  lie  in  a straight  line 
only  if  they  are  projected  down  onto  the  a-8  plane  in  direction  cosine 
space.  It  is  therefore  clear  that  varying  the  angle  of  incidence  merely 
shifts  the  diffracted  wave  field  in  direction  cosine  space  without 
changing  its  functional  form.  This  has  been  verified  experimentally  by 
mounting  a diffraction  grating  at  the  center  of  a transparent  hemisphere, 
placing  graph  paper  cn  the  plane  of  the  diffraction  grating  (a-8  plane), 
and  scribing  appropriate  latitude  lines  on  the  hemisphere  upon  which  the 
diffracted  ordt  . fall  when  illuminated  with  a small  laser  beam. 

Summary 

We  have  developed  a very  general  treatment  of  scalar  diffraction 
theory  that  yields  more  insight  and  accuracy  than  the  conventional 
Kirchhoff  and  Rayleigh- Sommerfeld  theories. 

By  describing  the  diffraction  process  in  terms  of  the  direction 
cosines  of  the  propagating  light  we  have  obtained  the  extremely  powerful 
result  that  the  diffracted  wave  field  on  an  observation  hemisphere  is 
given  directly  by  the  Fourier  transform  of  the  aperture  function.  This 
allows  us  to  apply  the  well-known  techniques  of  linear  systems  theory 
that  have  proven  so  useful  in  the  area  of  image  formation. 
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Furthermore,  we  have  shown  that  any  departures  of  the  actual 
diffracted  wave  field  from  that  predicted  by  tne  Fourier  transform  rela- 
tionship take  the  form  of  conventional  aberrations  whose  behavior  is 
well  understood  in  terms  of  the  dimensions  of  the  diffraction  aperture, 
the  radius  of  the  observation  hemisphere,  and  the  appropriate  field 
parameters. 
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CHAPTER  III 

SURFACE  SCATTER  THEORY 

In  the  following  treatment  the  scattering  of  light  from  optical 
elements  is  considered  to  be  solely  a surface  phenomenon.  Although  it 
is  recognized  that  bulk  scattering  mechanisms,  such  as  photon-phonon 
interactions  [Bloembergen  (1965)]  and  scattering  from  free  electrons 
[Vachaspati  (1964)],  can  exist  if  the  substrate  material  is  not  per- 
fectly conducting.  Leader  and  Dalton  (1972)  have  obtained  some  measured 
data  on  the  bulk  scattering  from  dielectrics.  The  excitation  of  surface 
plasmons  has  also  been  suggested  by  several  investigators  as  contribu- 
ting to  short  wavelength  scattering  from  polished  metal  surfaces 
[Beaglehole  (1970),  Hunderi  and  Beaglehole  (1970),  Crowell  and  Ritchie 
(1970),  Elson  and  Ritchie  (1971),  and  Daudd,  Savary,  and  Robin  (1972)]. 
However,  the  above  effects  are  believed  to  be  small  for  most  visible  and 
infrared  radiation  scattered  from  metal  surfaces. 

Surface  Scatter  Phenomena  as  a Diffraction  Process 
In  Chapter  II  it  was  shown  that  the  diffracted  wave  field  on  a 
hemisphere  is  given  directly  by  the  Fourier  transform  of  the  complex 
amplitude  transmittance  of  the  diffracting  aperture.  In  general,  the 
pupil  function  in  a diffracting  apert’ve  will  include  perturbations  that 
lead  to  scattered  radiation.  Surface  scatter  phenomena  can  thus  be 
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dom  phase  variations  in  addition  to  any  existing  amplitude  variations. 
The  diffraction  theory  of  the  previous  chapter  can  therefore  be  directly 
applied  to  the  problem  of  predicting  the  complex  amplitude  distribution 
on  an  observation  hemisphere  of  radius  r resulting  from  an  incident  beam 
of  light  being  reflected  from  a rough  surface. 
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Effect  of  the  Scattering  Surface  upon  the  System  Pupil  Function 

Figure  10  illustrates  the  surface  height  variations  (Wr)  as  a 
function  of  distance  along  the  surface.  This  surface  profile  has 
associated  with  it  an  autocovariance  function  and  a surface  height  dis- 
tribution function  as  indicated. 

Tne  effect  of  the  rough  surface  can  be  considered  to  be  a space- 
dependent  modifier,  or  random  component,  of  the  effective  pupil  function 
o*  the  system.  Hie  disturbance  emerging  from  the  scattering  surface  is 
given  by 


Uq(x  y;0)  = a^pCS.i/jO)  = aj?L(x,y;0')pR(3:,y;0) . (31) 

Here  the  pupil  function  of  the  system  producing  the  incident 
beam  is  given  by 

i2tiWL(x,y;Q ) 

PL&.y]  0)  = aL(x,y; 0)  e , (32) 

where  a ^ describes  the  amplitude  variations  across  the  exit  pupil  of  the 
system,  and  U ^ describes  any  phase  variations  or  aberrations  in  the 
wavefront  of  the  incident  beam. 
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Fig.  10.  Illustration  of  Surface  Heigi t Variations  and  Associated 
Statistical  Parameters. 


The  random  component  of  the  pupil  function  due  to  the  scattering 
surface  similarly  has  an  amplitude  and  phase  component 

i^vWR  (x,y;0) 

PH(x,y; 0)  = e , (32a) 

where  E is  the  reflectance  of  the  scatteiing  surface  and  is  the 
above-mentioned  surface  height  variation.  Note  that  the  phase  varia- 
tions on  the  reflected  wavefront  are  twice  as  large  as  the  actual 
variations  on  the  reflecting  surface. 
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Scattered  Wave  Field  on  a Hemisphere 

The  results  of  the  previous  chapter,  [Eq.  (29)],  indicate  that 
the  complex  amplitude  distribution  on  the  observation  hemisphere  is. 
given  directly  by  the  Fourier  transform  of  U0(x,y; 0) 

£2n& 

0(o,e;iO  = Y S-r-  Wo(a»S;0)}.  (33) 

iv 


The  total  reflected  flux  P is  obtained  by  applying  Rayleigh's 
theorem  from  Fourier  transform  theory 


00  00 
p = jj  |0(a,e;*)|zcM3  = a0z  Jj|p(x’,p';0) \2dx'dy'. 

-«®  (34) 

Noting  that  du>  = da,d&/y,  the  radiant  intensity  of  the  scattering 
system  can  be  written  as 

■T(a»6)  = = ~ iy(a,6;r)|2  = ao2y\&{p(.x,y;0) } |2.  (35) 

Utilizing  the  autocorrelation  theorem  of  Fourier  transform  theory,  this 
is  equivalent  to 


I(o,8) 


azy&{  P(£',p';0)p*(£'-£,p'-p;0)dr'£$'>.  (35a) 


Tne  autocorrelation  function  of  a Lambertian  surface  approaches 
a Dirac  6-function;  hence,  its  Fourier  transform  is  constant  and  we 
obtain 

■*  yX,  K - constant  (35b) 

which  is  consistent  with  Lambert's  cosine  law. 
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The  Effective  Transfer  Function  of  a Scattering  System 

Following  the  standard  procedure  used  in  image  evaluation,  the 
effective  transfer  function  of  the  scattering  system  is  defined  as  the 
normalized  autocorrelation  of  the  pupil  function 


o 


u 


0)  = 


co 

JJ  p(x',y'  ;0)p*(ffl'-s>y,-p;0)<&'<2^/, 

.00 

CO 

JJ  IpCife* ,p* ;0> 


(36) 


The  effective  spread  function  of  the  scattering  system  is  now 
defined  in  the  usual  way  as  the  Fourier  transform  of  the  transfer 
function 


S(a,B)  = &{Jf(x,y;  0)}  = 


rjjj  f(«,,p,;0)p*(^,-x,p'-^;0)c£b'd^,| 

x ~uu  / 


w 

JJ  |p(£’*j/’;0) \2dx'dy' 


(37) 


Direct  substitution  from  Eq.  (34)  and  Eq.  (35a)  results  in  the  following 
expression  for  the  effective  spread  function  in  terms  of  the  radiant 
intensity  of  the  scattering  system 


5(<*»2)  = ^p-I(a,e). 


(37a) 


Scattering  from  Optical  Surfaces 
Let  the  height  variations  of  a given  illuminated  area  on  a 
scattering  surface  be  a two-dimensional  sample  function  W^(x,y  ;0) . A 
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random  process,  W^(,x,y;Q),  is  made  up  of  an  ensemble  of  such  functions 
as  shown  in  Fig.  11.  For  fixed  spatial  coordinates,  f/j?(*i,pi;0)  is  a 
random  variable . And,  for  a specific  sample  function  with  fixed  spatial 
coordinates,  is  a single  number. 


"itf.y) 


Fig.  11.  An  Cnscnblc  of  Two-Dimensional  Sample  Functions  h'j(x.y) 

Representing  Surface  Height  Variations  Const it!itcs  a Random 
Process  iii(x,y) . 

Two  random  variables,  Sr(xs  ,yj)  and  W(x2,y2)  with  fixed 
spatial  coordinates  arc  also  shown. 
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Assumptions  Concerning  the  Statistical  Properties  of  an  Optical  Surface 

We  are  primarily  interested  in  the  scattering  behavior  of  opti- 
cal surfaces.  The  following  assumptions  are  made  concerning  the  statis- 
tical properties  of  an  optical  surface  prepared  by  conventional 
fabrication  techniques  on  ordinary  optical  materials: 

1.  The  reflectance  h is  constant  over  the  entire  surface.  This 
assumption  is  not  essential  lrut  it  is  reasonable  and  fur- 
thermore it  has  been  shown  by  Shack  (1967)  that  phase  fluc- 
tuations will  dominate  over  amplitude  fluctuations  in  their 
combined  effect  on  the  spread  function. 

2.  Wfl(x,y; 0)  is  a single-valued  gaussian  random  process. 

3.  0)  is  at  least  locally  stationary  in  the  statistical 
sense  (i.e.,  surface  is  homogeneous  and  isotropic). 

4.  The  random  variables  W^{xi,y\; 0)  and  Wpix2,yzl^) » produced 
by  any  two  fixed  pairs  of  spatial  coordinates,  are  jointly 
normal  regardless  of  their  correlation  coefficient. 

5.  ^(x,y;0)  is  weakly  ergodic  (i.e.,  the  mean  and  autocorrela- 
tion function  determined  by  space  averages  using  a singJe 
sample  function  WJ).{x,y;0')  are  the  same  as  those  determined 

■AX' 

by  ensemble  averages). 
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1116  Effective  Transfer  Function  of  a Scattering  Surface 
Substituting  Eq.  (31)  into  Eq.  (36)  we  obtain 


oo 

w 


A A 


12t\[Wii-Wl2\  ■MnUfR2-wR2] 
aLlaL2  e e dx'dy' 


l< 

J . 


12  dx'dy ' 


(38) 


where 


aLI  = a^Cx'  ,y'  ;0) 


L2 


a^x’-x,  y'-y;  0) 


wli  * 

^2  = 0'-$;  °) 

WR1  = ^(x»,y*;  0) 

^2  = ^(x’-x,  fi’-y,  0), 


(39) 


The  above  expression  for  the  transfer  function  contains  the  ran- 
dom variables  WR^  and  WRR;  therefore,  taking  the  expected  value  we  have 


aLlaL2e 


r{JtT(x,y;0)}  = 


i2v[WL1-WL2\  ( U*[Wr2'-Nr2}) 

T I 


dx'dy  * 


(40) 


la£l2<^'dy ' 
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Since  the  random  variables  involved  are  assumed  to  be  stationary, 
the  expected  value  under  the  integral  is  independent  of  ic'  and  y'  and  y 
can  be  taken  outside  the  integral 


E{Jf(x,y;0)}  = E< 


iA-n[WRi-WR2] 


jj  aLlal 


i2n[WLl-^L2\ 


dx'dy' 


co 

jj  \aL\2dx'dy' 


We  now  recognize  the  normalized  autocorrelation  function  in  the  above 
equation  to  be  the  transfer  function  of  the  optical  system  producing  the 
incident  beam 


■f  i2v[Wir^L2^ 

')aLlaL2* 


jeL&,y;  0)  = 


w 

j \aL\2dx'dy' 


The  average  quantity  in  Eq.  (41)  is  therefore  the  equivalent  transfer 
function  of  the  scattering  surface 


Jtf  (x,£;0)  = £<e 


WR1-VR2] 


But  this  is  merely  the  joint  characteristic  function  (Papoulis 
(1965),  p.  225)  of  the  two  random  variables  W^j  and  Since  W ^ and 

W n are  jointly  normal  random  variables,  it  can  be  shown 

t\u 
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that  (see  Papoulis  (1965),  p.  226)) 


( ^4tt IWR2-WR2] 

£je 


= ei4ir(ni-n2)e-8iT2(o12-2C12+022) 


(44) 


where 


C12  = S{(^i-ni)(^2-n2)> 


(45) 


is  the  covariance  function  of  the  random  variables  and  Wbo.  But 

HI  Hu  HI 

and  are  identical  functions  merely  displaced  from  one  another;  hence 


oi  - o2 


aW' 


ni  = n 2 


(46) 


and 


^12  = £j^(£»y) 


autocovariance  of  Wn. 

H 


(47) 


The  equivalent  transfer  function  of  the  scattering  surface  is  thus  given 
by 


.#^(x,y;0) 


-(4to#)2 


1 


e 


Cfr&>y) 


(48) 


where  o^2  is  the  variance  of  the  surface  height  distribution  function 
and  C^(x,y)  is  the  two-dimensional  autocovariance  function  of  the 
surface. 

Considerable  insight  into  the  scattering  process  can  now  be 
obtained  by  considering  the  nature  of  this  transfer  function.  The  auto- 
covariance function  approaches  the  value  o^2  as  the  displacement 
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approaches  zero.  The  equivalent  transfer  function  thus  approaches  unity 


as  expected.  As  the  displacement  approaches  infinity  the  auto covariance 


function  approaches  zero  and  the  equivalent  transfer  function 


approaches  a plateau  of  height  exp[-(2irc£)2] . 


The  equivalent  transfer  function  of  the  scattering  surface  can 


thus  be  regarded  as  the  sum  of  a constant  component  and  a bell -shaped 


component  as  shown  in  Fig.  12a.  Equation  (48)  can  therefore  be  rewrit- 


ten as 


•#R{x,y-,  0)  = A + BQ(2t,y;  0), 


where 


A = e 


(4irc^): 


B - 1 - e 


-(4tto£0: 


0 

vr? 


£(x>y;0)  = 


(4tt  Of,)2 
e - 1 


The  Effective  Spread  Function  of  a Scattering  Surface 


The  significance  of  this  interpretation  of  the  equivalent  trans- 


fer function  of  the  scattering  surface  is  dramatically  shown  by  the 


inferred  properties  of  the  corresponding  spread  function.  Since  the 
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transfer  function  is  the  sum  of  two  separate  components,  the  equivalent 
spread  function  of  the  scattering  surface  is  the  sum  of  the  inverse 
Fourier  transforms  of  the  two  component  functions.  Tho  constant  compo- 
nent transforms  into  a delta  function,  and  the  bell-shaped  component 
transforms  into  another  bell-shaped  function  as  shown  in  Fig.  12b. 

Hence  the  scattering  surface  reflects  an  incident  beam  of  light  as  a 
specularly-reflected  beam  of  diminished  intensity  surrounded  by  a halo 
of  scattered  light.  Furthermore,  the  relative  power  distribution 
between  the  specular  component  and  the  scattered  component  of  the 
effective  spread  function  are  given  by  the  quantities  A and  B 
respectively. 

Note  that  as  more  and  more  light  is  scattered,  energy  is  trans- 
ferred from  the  specular  component  of  the  spread  function  into  the 
scattered  component  of  the  spread  function.  For  a perfectly  Lambertian 
reflector  the  specular  component  disappears  completely  from  the  spread 
function. 

Relationship  between  the  Spread  Function  and  the  BRDF 

From  the  above  analysis  it  is  clear  that  the  Bidirectional 
Reflectance  Distribution  Function  (BRDF),  (see  Appendix  C) , which  is 
appearing  frequently  in  the  radiometry  literature,  is  merely  an  infinite 
family  of  two-dimensional  spread  functions. 

In  general,  the  above  scattering  function  (or  the  BRDF  of  dif- 
fuse reflectance  measurements)  will  change  with  angle  of  incidence  just 
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as  the  point  spread  function  of  an  optical  imaging  system  will,  in 
general,  vary  with  the  field  position  of  the  point  source.  However,  the 
analysis  of  imaging  systems  is  greatly  simplified  by  assuming  an  iso- 
planatic  system  in  which  the  point  spread  function  does  not  change  with 
field  position  (and  this  is  a reasonable  assumption  for  many  practical 
imaging  systems).  Similarly,  the  analysis  of  light  scattering  systems 
will  be  greatly  simplified  if  they  can  be  shown  to  be  shift-invariant 
(i.e.,  if  the  shape  of  the  scattering  function  does  not  change  with  the 
angle  of  incidence).  If  this  is  true,  the  four-dimensional  BRDF  degen- 
erates into  a single  two-dimensional  spread  function. 

Shift  Invariance  in  Direction  Cosine  Space 

The  scattered  light  distribution  on  an  observation  hemisphere 
will  appear  to  consist  of  the  sum  of  two  components,  a core  which  is  the 
delta  function  convolved  with  the  spread  function  of  the  optical  system 
producing  the  incident  beam,  and  a scattering  function  which  is  the 
bell-shaped  halo  convolved  with  the  spread  function  of  the  optical 
system. 

In  Fig.  13  we  have  merely  replaced  the  diffracting  aperture  of 
Fig.  8 with  a scattering  surface  and  the  geometry  of  the  measurements 
has  been  folded  about  the  reflecting  plane.  Hence,  we  have  the  incident 
beam  striking  the  scattering  surface  at  seme  angle  of  incidence,  a 
specularly-reflected  beam  striking  the  observation  hemisphere,  and  the 
scattered  light  distribution  being  sampled  at  an  arbitrary  point  with 
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g.  13.  Geometrical  Relationship  Between  Incident  Beam,  Scatter 
ing  Surface,  the  Measured  Scattered  Light  Distribution, 
and  the  Resulting  Spread  Function. 


direction  cosine  coordinates  a and  3.  The  scattered  light  distribution 
on  the  hemisphere  will,  in  general,  change  shape  drastically  with  angle 
of  incidence— becoming  quite  skewed  and  asymmetrical  at  large  angles  of 
incidence.  However,  our  theory  predicts  that,  for  certain  surfaces  with 
well-behaved  statistics,  if  the  data  collected  on  the  hemisphere  is 
plotted  as  a function  of  the  direction  cosines  of  the  position  vector  of 
the  observation  point,  this  new  scattering  function  will  not  change 
shape  but  will  merely  be  shifted  in  direction  cosine  space  with  changes 
in  angle  of  incidence.  This  is  a rather  significant  development  which 
has  profound  implications  regarding  the  quantity  of  data  required  to 
completely  characterize  a scattering  surface.  However,  it  remains  to  be 
shown  whether  scattering  surfaces  of  practical  interest  can  experimen- 
tally be  shown  to  obey  these  predictions. 
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CHAPTER  IV 

SURFACE  SCATTER  MEASUREMENTS 
Apparatus 

An  instrument  has  been  designed  and  built  at  the  Optical  Sci- 
ences Center  for  making  scattered  light  measurements  op  a hemisphere  as 
described  in  the  previous  chapter.  A schematic  diagram  of  this  appara- 
tus is  shown  in  Fig.  14.  The  incident  light  passes  through  a chopper  so 
that  synchronous  detection  with  a PAR  lock-in  amplifier  can  be  made. 

The  mechanical  apparatus  shown  in  Fig.  15  is  located  in  a small  photome- 
tric darkroom  in  which  the  experiment  is  conducted.  A movable  arm  with 
folding  mirrors  can  be  positioned  to  direct  the  incident  beam  onto  the 
sample  at  any  desired  angle.  A lens  positioned  on  this  arm  focuses  the 
incident  radiation  onto  the  hemisphere  mapped  out  by  the  detector;  hence, 
the  geometrical  configuration  is  consistent  with  that  illustrated  in 
Fig.  13  of  the  previous  chapter.  Two  separate  driving  mechanisms  allow 
us  to  measure  the  scattered  light  distribution  over  the  entire  hemi- 
sphere bounded  by  the  plane  of  the  sample. 

The  Light  Source 

The  light  source  employed  is  a Spectra-Physics  Model  165  Argon 
Ion  Laser.  The  laser  is  operated  with  a light-regulated,  single- 
frequency output  which  assures  intensity  regulation  to  within 
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Fig.  15.  Apparatus  for  Measuring  Scattered  Light  Distribution  from 
Optical  Surfaces. 

(1)  Precision  rotary  table,  (2)  worm  gear  drive  for  arm 
supporting  detector,  (3)  photomultiplier  tube  and  fiber- 
optic probe,  (4)  sample  holder,  (5)  movable  arm  with  fold- 
ing mirrors  and  lens  for  directing  and  focusing  incident 
beam,  (6)  P.A.R.  lock-in  amplifier,  (7)  high  voltage  power 
supply  for  PMT. 
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one  percent.  The  measurements  were  made  with  approximately  20  mw  of 
power  in  the  incident  beam  at  a wavelength  of  0.5145  urn. 

The  Detector  Unit 

The  detector  is  a Phillips  one- inch,  end-on  photomultiplier 
having  an  S-20  photocathode.  Light  reaches  the  photomultiplier  by  way 
of  a rigid  fiber-optic  probe.  Such  a probe  offers  several  distinct 
advantages  in  light  sampling.  In  addition  to  allowing  increased  angular 
resolution  throughout  the  sampling  space,  and  enabling  us  to  sample 
within  half  a degree  of  the  incident  or  specularly- reflected  beams,  it 
provides  the  ability  to  control  the  field  of  view  of  the  detector  for 
the  purpose  of  stray  light  rejection. 

The  original  configuration  consisted  of  a rigid  fiber-optic 
bundle  bent  such  that  one  end  was  pointed  toward  the  illuminated  spot  on 
the  sample.  This  end  of  the  bundle  thus  acted  as  the  collecting  aper- 
ture for  the  detection  system.  The  other  end  of  the  fiber-optic  bundle 
protruded  into  the  photomultiplier  tube  housing  followed  by  a series  of 
baffles  to  limit  the  field  of  view  of  the  detector  as  shown  in  Fig.  16a. 
This  resulted  in  a detector  response  which  had  a gauss i an  dependence 
upon  field  angle. 

However,  by  introducing  a small  collecting  lens  and  a field  stop 
in  front  of  the  fiber-optic  bundle  as  shown  in  Fig.  16b,  the  baffles  can 
be  eliminated  and  a well-defined  field  of  view  of  any  desired  size  can 
be  obtained  by  properly  choosing  the  size  of  the  field  stop.  This  is 
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Fig.  16.  (a)  Previous  Configuration  of  Detector  Probe  Unit, 

(b)  New  Configuration  of  Detector  Probe  Unit. 


more  clearly  illustrated  in  Fig.  17.  A coated  doublet  with  a \0-mm 
focal  length  was  edged  down  and  mounted  in  a black  anodized  brass  tube 
3 mm  in  diameter.  A field  stop  allowing  a 5-degree  field  of  view  was 
fabricated  and  inserted  into  the  tube  at  the  rear  focal  plane  of  the 
lens.  This  assembly  was  then  positioned  onto  the  end  of  the  fiber-optic 
bundle.  The  detector  response  from  a small  (point  source)  light  source 
was  then  recorded  as  a function  of  field  angle  for  both  detector  probe 
configurations.  The  results  are  displayed  for  comparison  in  Fig.  18. 
Both  the  flat  response  and  the  sharp  cutoff  obtained  with  the  modified 
unit  are  highly  desirable  features.  The  flat  response  promises  to  elim- 
inate signal  variations  due  to  slight  misalignment  or  wobble  in  the 
mechanical  instrument  while  scanning  over  the  hemisphere.  The  ability 
to  keep  the  field  of  view  small  with  a very  sharp  cutoff  is  essential 
for  stray  light  rejection. 

The  scattered  light  flux  from  a polished  surface  varies  by  sev- 
eral orders  of  magnitude  over  the  angular  range  to  be  measured.  Hence 
the  linearity  of  the  PMT  was  measured  using  a calibrated  neutral  density 
wedge  and  several  known  neutral  density  filters  to  vary  the  incident 
flux.  The  resulting  linearity  curve  is  shown  in  Fig.  19  and  indicates  a 
deviation  of  less  than  1%  over  a range  of  five  orders  of  magnitude  of 
the  incident  flux. 
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The  Scanning  Mechanism 

The  mechanical  apparatus  for  measuring  the  angular  distribution 
of  light  scattered  from  optical  surfaces  was  shown  in  Fig.  15.  The 
detector  probe  unit  is  mounted  on  a rigid  arm  that  can  be  rotated  ir. 
either  of  two  orthogonal  directions.  Rotation  about  a vertical  axis  is 
accomplished  by  means  of  a massive  precision  rotary  table.  The  rigid 
arm  is  attached  to  the  rotary  table  by  means  of  a worm  gear  arrangement 
that  allows  rotation  about  the  horizontal  axis.  The  intersection  of 
these  two  axes  defines  the  sample  location. 

The  Incident  Beam 

A large  movable  arm  with  appropriate  folding  mirrors  is  used  to 
direct  the  incident  beam  onto  the  sample  at  any  desired  angle  of  inci- 
dence. Immediately  following  thu  last  folding  mirror  the  beam  is 
focused  onto  a pinhole  which  acts  as  a spatial  filter  as  shown  in 
Fig.  20.  Lens  L2  then  forms  an  image  of  the  pinhole  upon  the  hemisphere 
mapped  out  by  the  collecting  aperture  of  the  scanning  fiber-optic  probe. 
This  spatial  filter  assembly  eliminates  from  the  beam  incident  upon  the 
sample  any  light  scattered  from  the  folding  mirrors  as  well  as  any  dif- 
fraction effects  from  the  chopper  blade. 

Experimental  Procedures 

The  instrument  described  in  the  previous  section  allows  us  to 
measure  the  scattered  light  distribution  over  the  entire  hemisphere 
bounded  by  the  plane  of  the  sample  for  all  angles  of  incidence.  However, 
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in  order  to  limit  the  quantity  of  data  to  be  collected,  the  scattered 
radiation  field  will  be  sampled  in  two  principal  directions.  These 
include  the  entire  plane  of  incidence  and  a plane  containing  the  inter- 
section of  the  specular  beam  with  the  unit  hemisphere,  which  is  perpen- 
dicular to  both  the  plane  of  incidence  and  the  plane  of  the  sample 
being  tested  (see  Fig.  21).  This  particular  sampling  procedure  was 
chosen  because  each  sampling  direction  then  involves  one  fixed  coordi- 
nate in  direction  cosine  space  as  discussed  in  the  theory  section  of 
this  report.  Furthermore,  the  apparatus  was  designed  such  that  each  of 
the  two  independent  drive  mechanisms  corresponds  directly  to  a given 
coordinate  in  direction  cosine  space.  Hence  for  a given  observation 
point  determined  by  the  angles  6 and  $ displayed  on  the  apparatus,  the 
corresponding  coordinates  in  direction  cosine  space  are  given  by 

a = cose  sin# 

8 = sin6.  (51) 

Measurements  will  be  made  to  within  one  degree  of  the  specular 
beam  in  both  directions.  These  measurements  will  be  made  for  a variety 
of  angles  of  incidence  on  a set  of  samples  consisting  of  an  assortment 
of  glass  types  having  been  polished  to  varying  degrees  of  surface 
roughness.  A sample  with  a diffuse  MgO  coating  will  be  used  as  a refer- 
ence. Asymmetries  that  arise  as  a result  of  polishing  techniques, 
cleaning,  and  coating  anomalies,  can  be  investigated  by  rotating  the 
sa<4p;e  about  a normal  to  the  surface. 
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plane 


21.  Geometrical  Configuration  of  Two  Principal 
Plants  in  which  the  Scattered  Light  Field 
Will  Be  Sampled. 


Sample  Preparation 

An  original  goal  of  the  research  reported  here  was  to  examine 
the  scattering  properties  of  samples  whose  surface  characteristics  span 
those  typically  produced  with  optical  fabrication  techniques. 

Surface  preparation  techniques  used  to  produce  a set  of  samples 
are  outlined  in  Table  2.  All  samples  were  finished  to  be  nominally  flat. 

The  prepared  samples  were  cleaned  prior  to  coating  with  aluminum. 
Cleaning  consisted  of  careful  washing  with  Liquinox,  a mild  detergent, 
under  very  warm,  filtered  tap  water.  Samples  were  then  mounted  in  a 
sample  holder  while  held  in  distilled  water.  Once  in  a holder,  samples 
were  moved  to  an  ultrasonic  cleaner  filled  with  distilled  water  for 
rinsing.  Once  rinsed  the  samples  were  set  to  dry  in  a dust- free  atmo- 
sphere. Dry  samples  were  removed  from  the  sample  holders  and  placed  in 
individual  boxes  being  supported  by  the  edge  of  their  backside.  Mott 
(1S71)  used  a similar  cleaning  technique,  which  he  describes  more 
completely. 

Cleaned  dry  samples  were  then  placed  in  a high  vacuum  chamber 
and  coated  to  i ear  opacity  with  pure  aluminum.  Coating  technique  varied 
from  the  standard  only  in  that  excessive  care  was  taken  to  allow  the 
chamber  to  reach  a pressure  below  2 * 10‘e  torr  prior  to  coating.  The 
samples  were  allowed  to  cool  to  room  temperature  prior  to  removal  from 
the  chamber.  Each  coating  run  contained  ten  different  samples.  Once 
coated,  samples  were  returned  to  their  individual  storage  boxes.  After 
all  samples  were  coated,  the  best  samples  of  each  type  were  selected  for 
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measurement.  This  selection  was  made  on  the  basis  of  individual  inspec- 
tion of  each  sample  while  held  under  a microscope  illuminator  i.  an 
otherwise  dark  room.  Samples  that  had  coating  nonuniformities,  sleeks 
or  pinholes  were  rejected  as  were  those  with  waterspots,  large  scratches, 
or  otherwise  questionable  appearance.  Prior  to  each  set  of  scatter  mea- 
surements, samples  were  again  individually  inspected  for  flaws.  Dust 
was  removed  using  a commercially  available  pressurized  air  can.  After 
the  scattering  characteristics  of  each  sample  were  measured,  the  sample 
was  returned  to  its  individual  box. 

System  Alignment 

Prior  to  making  any  meaningful  scatter  measurements  it  was 
necessary  to  systematically  align  the  entire  system. 

The  incident  laser  beam  was  first  made  accurately  horizontal. 
Then  the  mechanical  apparatus  was  positioned  such  that  the  axis  of  rota- 
tion of  the  movable  arm  supporting  the  folding  mirrors  (see  Fig.  20)  was 
colinear  with  the  incident  beam.  This  was  accomplished  by  means  of  four 
massive  leveling  screws  at  the  base  of  the  stand  supporting  the  entire 
apparatus . 

With  mirror  Ml  removed  the  laser  beam  passes  through  the  small 
hole  in  the  center  of  the  bearings  upon  which  the  movable  arm  rotates. 
This  assures  that  the  incident  beam  is  indeed  colinear  with  the  axis  of 
rotation  and  furthermore  allows  the  sample  holder,  which  is  mounted  on  a 
shaft  passing  through  the  axis  of  the  precision  rotary  table,  to  be 
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accurately  positioned  so  that  the  center  of  the  scattering  surface  lies 
at  the  intersection  of  these  two  perpendicular  axes. 

Mirror  Ml  is  then  put  into  place  and  adjusted  until  the  laser 
beam  is  centered  upon  mirror  M2.  Similarly,  M2  is  adjusted  until  the 
beam  is  centered  upon  M3. 

With  the  lenses  LI  and  L2  and  the  pinhole  P removed  from  the 
system,  mirror  M3  is  adjusted  until  the  beam  is  centered  upon  the  sample 
holder. 

The  incident  beam  is  then  made  accurately  perpendicular  to  the 
axis  of  rotation  of  the  movable  arm  by  the  following  procedure.  With 
the  sample  removed,  the  rotary  table  is  positioned  so  that  the  beam  is 
centered  on  the  detector  at  position  A.  A polished  sample  is  placed  in 
the  holder  and  adjusted  until  the  specular ly-relfected  beam  returns  pre- 
cisely along  the  incident  beam.  Now  the  table  is  rotated  until  the 
detector  is  centered  upon  the  reflected  beam  at  position  B.  If  A and  B 
are  not  precisely  180°  apart,  systematically  adjust  mirrors  M2  and  M3 
and  repeat  the  above  procedure. 

Lenses  LI  and  L2  are  then  placed  in  the  beam  and  properly  cen- 
tered. And  finally,  the  pinhole  P is  accurately  positioned  at  the  back 
focal  position  of  lens  LI. 

Kith  the  system  properly  aligned,  the  movable  arm  can  be  rotated 
to  direct  the  inciaent  beam  at  any  desired  angle  without  requiring  other 
adjustments  to  keep  the  beam  centered  upon  the  sample. 
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An  additional  requirement  is  that  the  PMT  with  its  associated 
fiber-optic  probe  be  positioned  and  aligned  such  that  the  field  of  view 
of  the  PMT  remain  accurately  centered  on  the  illuminated  portion  of  the 
sample  throughout  the  entire  range  of  its  scanning  motion.  Provisions 
were  therefore  made  to  allow  three  degrees  of  freedom  (one  translation 
and  two  rotation)  in  adjusting  the  position  and  orientation  of  the  PMT 
housing.  Removing  the  PMT  from  its  housing  and  illuminating  the  fiber- 
optic probe  from  the  back  side  greatly  facilitates  this  alignment  pro- 
cedure as  it  allows  one  to  directly  observe  the  field  of  view  on  the 
sample  holder  while  making  the  necessary  adjustments. 

Measurement  Technique 

The  sample  to  be  measured  is  placed  in  the  holder  and  the 
movable  arm  positioned  to  achieve  the  desired  angle  of  incidence.  A 
calibrated  attenuator  is  then  placed  in  the  incident  beam  and  the  detec- 
tor centered  on  the  specularly- reflected  beam.  The  collecting  aperture 
of  the  fiber-optic  probe  is  large  enough  to  collect  the  entire  specular 
beam;  hence,  the  output  signal,  of  the  PMT  in  this  position  is  pro- 
portional to  the  total  flux  in  the  specular  beam. 

The  detector  • s then  moved  a known  angular  distance  (approxi- 
mately one  degree)  from  the  specular  beam  and  the  attenuator  removed.  A 
profile  of  the  scattered  light  distribution  is  then  measured  by  scanning 
the  observation  hemisphere  with  the  fiber-optic  probe.  Approximately 
30  separate  readings  are  taken  at  different  angular  positions  between 
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the  specular  beam  and  the  plane  of  the  sample.  This  set  of  readings  is 
the  raw  data. 

The  sample  is  then  removed  and  the  incident  beam  allowed  to  pass 
unobstructed  through  the  sample  holder  and  into  a black  absorbing 
Rayleigh's  horn.  Background  measurements  are  then  made  along  the  same 
profile  as  above  and  subtracted  from  the  raw  data.  These  background 
measurements  were  found  to  be  completely  negligible  in  most  instances. 

Tne  data  now  describes  the  spread  function  of  the  scattering 
system,  which  is  made  up  of  the  spread  function  of  the  scattering  sur- 
face convolved  with  the  spread  function  of  the  optical  system  producing 
the  incident  beam.  These  are  shown  in  Fig.  22. 

The  spread  function  of  the  incident  beam  is  then  measured  by 
again  placing  the  calibrated  attenuator  in  the  incident  beam  and  cen- 
tering the  detector  on  the  direct  beam  passing  through  the  empty  sample 
holder.  Since  the  collecting  aperture  of  the  fiber-optic  probe  is  large 
enough  to  collect  most  of  the  incident  beam,  the  output  signal,  of 
the  I’MT  in  this  position  represents  a good  approximation  to  the  total 
flux  in  the  incident  beam.  The  detector  is  then  moved  a known  angular 
distance  (approximately  one  degree)  from  the  direct  beam  and  the  attenu- 
ator removed.  A profile  of  the  incident  beam  is  then  measured.  These 
readings  rapidly  diminish  to  zero  within  five  degrees  of  the  peak  value. 

Since  the  spread  function  of  the  incident  beam  (Fig.  22a)  is 
narrow  compared  to  the  scattering  function  of  the  surface  (Fig.  22b), 
the  scattei ing  portion  of  the  surface  spread  function  is  virtually 
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unchanged  by  the  convolution  operation  while  the  delta  function  compo- 
nent merely  replicates  the  beam  spread  function.  The  desired  scattering 
function  can  thus  be  obtained  by  subtracting  the  beam  spread  functi  n 
readings  from  the  raw  data. 


Data  and  Results 

Rather  extensive  scatter  measurements  have  been  made  on  two 
representative  surfaces.  One  is  a ground  gla-s  surface  which  is  a very 
diffuse  reflector.  The  other  is  an  optically  polished  surface  which  is 
a nice  specular  reflector.  Both  samples  were  coated  with  aluminum  prior 
to  making  the  measurements.  Four  separate  scattering  profiles  from  the 
specular  beam  to  the  plane  of  the  sample  (see  Fig.  21)  were  made  at 
several  angles  of  incidence. 

The  backscattering  profile  of  the  scattered  light  distribution 
for  these  two  samples  is  shown  in  Fig.  23  for  several  different  angles 
of  incidence.  Note  that  we  are  plotting  the  quantity  V/  along  the 
ordinate.  Since  V is  proportional  to  the  power  collected  per  unit  solid 
angle  subtended  by  the  fiber-optic  probe  and  VD  is  proportional  to  the 
total  power  in  the  incident  beam,  we  can  write 

V P/Aw  _ P/ (4yAwj  er>4>r) 

" Pay  ~ vo/&  ~ 1 


whe  •»  4 is  the  illuminated  area  on  the  sample. 

Hence  :he  quantity  plotted  along  the  ordinate  is  equivalent  to 
the  reflected  radiance  in  the  sampled  direction  divided  by  the  incident 
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irradiance,  which  is  precisely  the  way  the  BRDF  is  defined.  Our  scat- 
tering curves  are  therefore  one-dimensional  profiles  of  the  four- 
dimensional  BRDF.  Also,  in  accordance  with  our  theory  we  are  plotting 
this  function  versus  the  quantity  B-p0,  which  is  the  distance  of  r.he 
observation  point  from  the  specular  beam  in  direction  cosine  space . 

Both  coordinates  are  then  plotted  on  a logrithmic  scale.  For  the  pol- 
ished sample,  the  five  curves  with  the  incident  angle  varying  from  zero 
to  60°  coincide  almost  perfectly.  Hence,  it  is  apparent  that  the 
scattering  function  does  not  appreciably  change  with  the  angle  of 
incidence.  The  corresponding  curves  for  the  rough  sample  coincide  for  a 
substantial  range  of  angles  then  begin  to  depart  somewhat  at  the  large 
angles. 

The  four  separate  profiles  of  the  scattered  light  distribution 
from  the  same  two  samples  with  the  incident  beam  at  45°  are  shown  in 
Fig.  24.  Again  the  curves  for  the  polished  sample  coincide  almost  per- 
fectly, s-  gesting  a rotationally-symmetric  distribution  in  direction 
cosine  space.  Some  asymmetry  is  noted  in  the  scattered  light  distribu- 
tion from  the  rough  sample. 

The  data  on  Figs.  23  and  24  confirm  that  for  a certain  class  of 
surfaces  (in  which  optically  polished  glass  is  definitely  a member,  and 
ground  glass  can  perhaps  be  included  to  a lesser  extent),  the  scattering 
properties  are  indeed  shift  invar: ant,  and  can  be  completely  character- 
ized by  a single  set  of  measurements  at  a fixed  angle  of  incidence.1 


The  scattering  function  for  a variety  of  samples  with  a wide 
range  of  ims  surface  roughnesses  are  exhibited  in  Fig.  25.  Note  that 
the  MgO  surface,  which  is  somewhat  of  a diffuse  reflectance  standard, 
yields  a straight  horizontal  line  as  a Lambertain  reflector  should. 

Since  the  radiance  of  a Lambertian  source  is  given  by  [Hudson  (1969)] 

L - M-  (55) 

where  M is  the  total  eraittance  into  a hemisphere,  the  height  of  the  MgO 
scattering  curve  is  taken  to  be  1/n  and  all  other  scatter  measurements 
are  normalized  accordingly.  The  values  of  the  rms  surface  roughness  are 
merely  estimates. 

It  may  be  of  interest  to  compare  the  scatteiing  curves  resulting 
from  special  materials  or  unusual  fabrication  techniques  with  those  of 
more  conventional  optical  surfaces.  For  example,  the  result  of  scat- 
tering measurements  on  a polished  beryllium  sample  and  an  ionically- 
polished  fused  quartz  sample  are  shown  in  Fig.  26  along  with  some  curves 
from  conventional  optical  surfaces.  All  measurements  are  normalized  by 
the  reflectance  of  the  surface  so  that  the  beryllium  does  not  appear  to 
have  better  scattering  characteristics  due  to  its  lower  reflectance. 
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Fig.  26.  Comparison  of  Scatter  Profiles  from  Special  Materials 
and  Unusual  Fabrication  'iechniques  with  Those  of 
Conventional  Optical  Surfaces. 


CHAPTER  V 


SURFACE  STRUCTURE  MEASUREMENTS 

Three  traditional  techniques  for  obtaining  surface  structure 
information  involve  profilometry,  electron-microscopy,  and  FECO  inter- 
ferometry. Most  profilometers  provide  too  coarse  a measurement  for 
optical  surfaces.  The  electron  microscope  works  nicely  on  the  rough 
ground  glass  surfaces  but  fails  to  yield  sufficient  information  about 
the  smooth  polished  surfaces.  However,  since  it  requires  multiple 
reflections,  the  FECO  interferometer  works  well  on  smooth  surfaces  with 
a strong  specular  beam  but  does  not  yield  good  results  for  the  rough 
diffuse  surfaces.  The  latter  two  complementary  techniques  will  thus  be 
utilized  in  our  research  effort. 

Electron  Microscopy 

Surface  profiles  of  the  rough  samples  are  determined  from 
electron-micrograph  stereo  pairs  using  conventional  stereo- 
photogrammetric  techniques  [Moffitt  (1959)].  Nankivell  (1963)  discusses 
some  of  the  stereo-photogrammetric  problems  unique  to  electron  micro- 
scope applications.  A typical  electron-micrograph  stereo  pair  is  shown 
in  Fig.  27a  with  a line  scribed  to  indicate  the  position  of  a set  of 
preliminary  surface  height  measurements  that  were  made  with  a standard 
Fairchild  Stereocomparagraph.  This  instrument  consists  of  a mirror 
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stereoscope  fitted  with  a parallax  bar  containing  a micrometer  for  mea- 
suring the  parallax  of  each  desired  pair  of  points.  The  resulting  sur- 
face profile  is  shown  in  Fig.  27b. 

The  two  statistical  parameters  which  determine  the  scattered 
light  characteristics  are  the  variance  of  the  surface  height  distribu- 
tion and  the  surface  autocovariance  function.  A computer  program  was 
written  which  takes  the  surface  profile  data  and  determines  the  above 
two  parameters.  The  transfer  function  of  the  scattering  surface  is  then 
calculated  from  £q.  (40)  of  Chapter  III.  The  Fourier  transform  of  this 
function  then  yields  the  predicted  scattered  light  distribution. 

A Hewlett-Packard  Model  9810  programmable  desktop  calculator  and 
plotter  is  used  for  this  data  reduction.  The  program  has  only  recently 
been  debugged  and  become  operational. 

Figure  28  shows  several  sections  of  a suvface  profile,  a histo- 
gram indicating  the  surface  height  distribution  along  with  a gaussian 
function  for  comparison,  the  surface  autocovariance  function,  the 
resulting  transfer  function,  and  finally  the  scattering  function  or 
spread  function. 

As  this  is  a very  recent  result  and  due  to  some  scaling  problems 
in  the  computer  program,  a comparison  of  this  predicted  spread  function 
with  the  measured  scattering  function  is  not  yet  available.  This  work 
is  being  continued  under  Contract  F94701-75-C-0106. 
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FECO  Interferometry 


Surface  profiles  of  the  smooth  samples  will  be  made  with  a FECO 
ir  erferometer.  This  data  will  also  be  used  to  determine  the  variance 
of  the  surface  height  distribution  and  the  surface  autocovariance  func- 
tion. Bennett  (1974)  has  reported  on  a scanning  FECO  interferometer 
that  can  measure  very  small  height  differences  with  a lateral  resolution 
of  2 microns  to  yield  statistics  for  optical  surfaces.  This  instrument, 
along  with  auxiliary  equipment  which  includes  a slow-scan  TV  camera, 
signal  averager,  mini-computer,  and  a teletype  unit,  yields  the  surface 
profile,  rras  roughness,  surface  height  distribution  function,  autocovar- 
iance function,  and  other  statistical  parameters  for  the  surface. 

We  have  sent  two  samples  to  be  measured  with  the  above  apparatus. 
Additional  samples  will  be  measured  on  a Hilger  Watts  N130  FECO  inter- 
ferometer here  at  the  Optical  Sciences  Center.  The  results  of  these 
measurements  will  be  used  to  predict  a surface  spread  function  by 
applying  the  theory  of  Chapter  III.  This  predicted  spread  function  will 
then  be  compared  to  the  directly  measured  scattering  function  in  an 
attempt  to  experimentally  verify  our  scattering  theory.  The  results  of 
this  work  will  be  reported  under  Contract  F04 701-75-C-0106. 


CHAPTER  VI 
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SUMMARY 

A theoretical  model  to  explain  surface  scatter  has  been  devel- 
oped in  which  the  scattering  process  is  assuaed  to  be  a diffraction 
phenomena. 

By  describing  the  diffraction  process  in  terms  of  the  direction 
cosines  of  the  propagating  light  we  have  obtained  the  extremely  powerful 
. jult  that  the  diffracted  wave  field  on  an  observation  hemisphere  is 
given  directly  by  the  Fourier  transform  of  the  aperture  function.  This 
allows  us  to  apply  the  well-known  techniques  of  linear  systems  theory 
that  have  proven  so  useful  in  the  area  of  image  formation.  Furthermore, 
we  have  shown  that  any  departures  of  the  actual  diffracted  wave  field 
from  tha.  predicted  by  the  Fourier  transform  relationship  take  the  form 
of  conventional  aberrations  whose  behavior  is  well  understood  in  terms 
of  the  dimensions  of  the  diffraction  aperture,  the  radius  of  the  obser- 
vation hemisphere,  and  the  appropriate  field  parameters. 

This  diffraction  theory  is  then  generalized  to  include  the  scat- 
tering effects  of  rough  surfaces.  For  certain  surfaces  with  well- 
behaved  statistics,  our  theory  predicts  that,  if  the  data  collected  on  a 
hemisphere  is  plotted  as  a function  of  the  direction  cosines  of  the 
observation  point,  this  new  scattering  function,  which  depends  only  upon 
the  variance  of  the  surface  height  distribution  and  the  surface 
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autocovariance  function,  will  not  change  shape  but  will  merely  be 
shifted  in  direction  cosine  space  with  changes  in  angle  of  incidence. 
This  is  a rather  significant  development,  which  has  profound  implica- 
tions regarding  the  quantity  of  data  required  to  completely  characterise 
a scattering  surface.  No  approximations  are  made  in  the  theoretical 
development  concerning  the  size  of  the  surface  height  variations;  hence, 
this  theory  should  work  equally  well  for  very  smooth,  highly-polished 
optical  surfaces  and  the  very  rough  surfaces  frequently  used  for  baffle 
materials. 

An  instrument  has  been  designed  and  built  at  the  Optical  Sci- 
ences Center  for  making  scattered  light  measurements  (or  BRDF  measure- 
ments) in  an  attempt  to  verify  the  above  theory.  This  instrument  is 
capable  of  making  measurements  over  the  entire  hemisphere  bounded  by  the 
plane  of  the  sample  (with  an  angular  resolution  of  less  than  one  degree! 
for  desired  angle  of  incidence. 

Extensive  scatter  measurements  at  a variety  of  angles  of  inci- 
dence have  been  made  on  two  representative  surfaces.  One  is  a ground 
glass  surface  that  is  a very  diffuse  reflector.  The  other  is  an  opti- 
cally-polished surface  that  is  a nice  specul.tr  reflector.  Both  samples 
were  aluminized  prior  to  making  the  measurements.  This  data  was  then 
plotted  as  a function  of  the  direction  cosines  of  the  position  vector  of 
the  observation  point.  For  the  case  of  the  smooth  sample,  the  five 
curves  with  the  angle  of  incidence  varying  from  zero  to  ou  degrees  coin- 
cide almost  perfectly.  The  scattering  curves  for  the  rough  sample 


coincide  for  a substantial  range  of  angles  then  begin  to  depart  somewhat 
at  large  ancles.  The  data  from  these  two  surfaces  confirm  that  for  a 
certain  class  of  surfaces  (in  which  optically-polished  glass  is  defi- 
nitely a member,  and  ground  glajs  can  perhaps  be  included  to  a lesser 
degree)  the  scattering  properties  of  the  surface  are  indeed  shift 
invariant,  and  can  be  completely  characterized  b£  £ single  set  of 
measurements  at  a fixed  ang?e  of  incidence. 

Surface  profile  measurements  have  been  made  on  the  rough  sample 
from  electron-micrograph  stereo  pairs  using  conventional  stereo- 
photogrammetric  techniques.  Similar  data  wilJ  be  obtained  from  a FECO 
interferometer  for  the  smooth  sample.  A computer  program  has  teen 
written  that  takes  the  surface  profile  data  and  determines  the  variance 
of  the  surface  height  distribution  and  the  surface  autocovariance  func- 
tion. These  two  statistical  parameters  are  then  used  to  calculate  the 
effective  transfer  function  of  the  scattering  surface.  The  Fourier 
transform  of  this  quantity  yields  the  predicted  spread  function  or  scat- 
tering function,  which  can  then  be  directly  compared  to  the  measured 
scattered  light  distribution.  Tnis  wot'.,  is  being  continued  under  Con- 
tract F04701-75-C-0106. 

The  above  developments  indicate  that  the  inverse  scattering 
problem  (determining  surface  properties  from  scattered  light  measure- 
ments) may  become  far  more  attractive  than  measuring  surface  properties 
directly. 
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APPENDIX  A 


FRESNEL  AND  FRAUNHOFER  DIFFRACTION 


In  order  to  simplify  the  calculation  of  diffracted  wave  fields, 
it  is  customary  to  impose  the  restriction  that  the  observation  distance 
z in  Eq.  (15)  be  much  larger  than  the  maximum  linear  dimension  of 
the  diffracting  aperture  or  the  region  of  interest  on  the  observation 
plane.  Under  this  condition  the  term  with  (1  + 6)2  in  the  denominator 
will  not  differ  significantly  from  unity  and  the  quantity  {Z  - 2)  in  the 
exponent  is  adequately  approximated  by  retaining  only  the  first  two 
terms  of  the  binomial  expansion  for  Z, 


Z - £ 


fx2+D*)  + (x,2+f/,2l 


(xx 1 +yy 1 ) 


(A.l) 


A sufficient  condition  for  the  validity  of  the  above  expression  is  the 
following  Fresnel  approximation 


s3  >>  7 [(x-x')2  + ty-y')2Y 


(A.  2) 


Substituting  Eq.  (A.l)  into  Eq.  (15)  results  in  the  familiar  Fresnel 
diffraction  integral 


86 


U(x,y;z ) = 


i2us  ^ (x2+y2) 


CO 
• * 

U0&',y';i 

J 


Gr,2+p'2) 


(xx'+yp') 


dx'dy' , 


(A.  3) 


which  is  valid  in  the  Fresnel  region  defined  by  Eq.  (A. 2). 

If  the  following  more-stringent  condition,  called  the 
Fraunhofer  approximation,  is  imposed  upon  the  observation  distance. 


2 >■•  ttC®'2  + y'2)max  , 


(A.  4) 


then  the  quadratic  phase  function  in  the  aoove  equation  is  approximately 
unity  over  the  entire  aperture.  The  diffracted  wave  field  on  an  obser- 
vation plane  in  the  Fraunhofer  region  defined  by  Eq.  (A. 4)  is  thus  given 
directly  from  the  Fourier  transform  of  the  pupil  function  U,-{x'  ,y' 0). 
This  is  the  well-known  Fraunhofer  diffraction  integral 


£2it3  (£2+y2)  f- 

U(x,y;z)  = — — e j U0(x'  ,y'  ;D) 

T2  ** 


2 7F  /A  A I A A 

- (XX'  + yy ') 

x e dx'dy'. 


(A.  5) 


However,  the  pupil  function  U0(x' ,y' ;0)  is  equal  to  the  product 
of  the  complex  amplitude  transmittance  of  the  diffracting  aperture 
(aperture  function)  T0(x' ,y ' ;0) , and  the  comolex  amplitude  distribution 
incident  upon  the  aperture.  Hence,  only  who-'  a plane  wave  is  incident 
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upon  the  aperture  is  the  Fraunhofer  diffraction  pattern  given  directly 
by  the  Fourier  transform  of  the  aperture  function 


U(x,y;z) 


i2nz  ij-  C x2+y 2) 
- e 3 

• A 

tz 


- ^ Cxxi+|/P') 

x s dib'dy' . (A. 6) 

Note  that  an  identical  expression  is  obtained  from  the  Fresnel  diffrac- 
tion formula  if  a spherical  wave  converging  to  the  observation  plane  is 
incident  upon  the  diffracting  aperture 

_ H.  (£»2+p»2) 

yo(x’,p’;0)  = T0@',y'; 0}  e * . (A. 7) 


In  both  of  the  above  cases  the  Fourier  transform  must  be  evaluated  at 
a = x/z  and  0 = x/z  to  assure  proper  space  scaling  in  the  observation 
plane. 
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APPENDIX  B 


CALCULATION  OF  ABERRATION  COEFFICIENTS 


FOR  DIFFRACTED  WAVE  FIELDS 


For  the  case  of  a plane  wave  incident  upon  an  aperture,  the 
diffracted  wave  field  on  an  observation  plane  is  given  by  Eqs.  (16)  and 
(17),  where 


W = (t-£)  + (x'x  + y’y)/z. 

The  quantity  Z can  be  written  as 


l = J tx-x  • ) 2 + 'T2  + 


v*  ,ere 


= z /I  + [§2+§»2  - 2(;ttfc,+0?Y)]/*2, 


a2  = x2  + y2,  s'2  ■«  x’2  + y'2. 


A binomial  expansion  of  the  above  square  root  results  in  the  following 
expression  for  W 


W = y [S*  + t>,2]/z2 


- | (s'*  + s’**  + 4(xx,+yy’}2  + 2 a2s'2  - 4s2(xx'+yy*) 


- 4s,2(xx'+yy,)]/sl4 


+ higher-order  terms. 
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Equating  coefficients  of  corresponding  terms  between  this  equa- 
tion and  the  wavefront  aberration  function  given  by  Eq.  (24),  we  obtain 
the  aberration  coefficients  tabulated  in  the  first  column  cf  Table  1, 
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If  we  now  have  a spherical  wave  incident  upon  the  aperture,  the 
quantity  W in  Eq.  (17)  is  given  by 


W = ( l-z ) - ( l0-z ) + (&&'  + 


(B.6) 


where 


;s 


& 


/r'2-  + p'2  + 22  = z /I  + 3'2/jg2. 


A binomial  expansion  of  this  quantity  results  in 


£0  - 2 = f (s'/z)2  - g (s’/2)4  + 


(B.7) 


But  these  will  merely  cancel  identical  terms  in  the  previous  expansion 
thus  resulting  in  the  aberration  coefficients  presented  in  column  two  of 


Table  1. 


For  the  c-se  of  a hemispherical  observation  space  the  diffracted 
wave  field  is  given  by  Eqs.  (22)  and  (23).  With  a plane  wave  incident 
upon  the  aperture,  we  have 


W - (£-r)  + (ax'  + 3p'). 


(B.8) 


The  quantity  t can  be  written  as 


l - /(£-*')2  + (£-£') 2 7 z2 


= f /T7~[§’?.  - 2Ho£'  + 8#’)]/f*2  , 


(B.9) 


where 


r'2  x*2  + y'2,  r2  = x2  + y2  + s2. 


A binomial  expansion  of  the  above  square  root  yields 


W 


= j is'/r)2  - ^ [s'4  - 4rs,2(c 'x'+Bp')  + 4r2(ax'+Bi/,)2]/r4 


+ higher-order  terms.  (B.10) 

If  we  again  assume  a rotationally-symmetric  diffracting  aperture 
we  can,  without  loss  of  generality,  choose  the  observation  point  on  the 
p-axis.  Let  us  therefore  set  a = 0.  We  can  also  let  y'  = s'  cos<|>  which 
results  in 


W = j (a'/^)2  - 7T  [s'4  - 4rBs'3  cos<J>  + 4r2B2s’2  cos2<J>3/r4 


+ higher-order  terms. 


(B. 11) 


If  we  now  substitute 


6 = PB. 


max’ 


A U 

S = a2 


into  the  previous  equation,  we  obtain 

& ~ § id/lr)2  az  - j [id/ 2) 4 a4  - 4rBmax  id/ 2) 3 Bp3  cos4> 

+ 4^%ax2  C<-*/2) 2 32p2  cos24>]/r4 

+ higher-oi'der  terms.  (B.12) 

Again  equating  coefficients  of  corresponding  terms  between  this 
equation  and  the  wr.^e front  aberration  function  given  by  Eq.  (24),  we 
obtain  the  aberration  coefficients  tabulated  in  column  three  of  Table  1. 
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If  we  now  have  a spherical  wave  illuminating  the  aperture  and  a 
hemispherical  observation  space,  the  quantity  V in  Eq.  (23)  is  given  by 

W = (1-r)  - ( l0-r ) + (otic*  + By'),  (B.13) 

where 

lQ  = ^2  + y'7'WW  = f*  A + §'2/f>2. 

A binomial  expansion  of  this  quantity  results  in 

lQ  - r = | (s«/r)2  - | (s'/r)4  + ...  . (3.14) 

Once  again  these  terms  merely  cancel  identical  terms  in  the  previous 
expansion,  leaving  only  coma  and  astigmatism  present  in  the  diffracted 
wave  field  as  indicated  in  the  last  column  of  Table  1. 
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APPENDIX  C 

BIDIRECTIONAL  REFLECTANCE  DISTRIBUTION  FUNCTION 

The  basic  quantity  that  characterizes  (geometrically)  the 
reflecting  properties  of  a surface  element  <34  is  the  bidirectional 
reflectance  distribution  function  (BRDF).  This  quantity 

= 6r, <PP;  (sr*1) 

(C-l) 

is  defined  by  Nicodemus  (19701,  as  the  reflected  radiance 

9r,^ fi'i)  of  the  surface  element  dA  in  the  direction  (0^,, 
divided  by  the  incident  irradiance  d&£(0£,<j>£)  = Li  producing 

it.  The  geometry  of  this  situstion  is  illustrated  in  Fig.  29,  where  the 
element  of  projected  solid  angle  is  given  by  d£l  = cosedw. 

Tin;  umerical  value  of  the  BRDF  for  a given  pair  of  incident  and 
reflected  ray  directions  may  vary  from  zero  to  infinity.  In  particular, 
consider  two  ideal  cases.  The  BRDF  is  a constant  *or  all  reflected 
directions  for  a perfectly  diffuse  (Lambertian)  surface;  and  it  becomes 
infinite  (as  a Dirac  delta  function)  for  a perfectly  specular  reflector. 
The  BRDF,  defined  above  as  a ratio  of  infinitesimals,  is  an  idealized 
concept  that  can  never  be  measured  exactly.  Real  measurements  are 
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always  made  over  some  finite  solid  angle  and  wavelength  interval  and  can 
therefore  yield  only  average  values  f?  over  those  parameter  intervals. 

The  BRDF  is  basic  in  the  sense  that  all  other  reflectance  or 
scattering  functions  can  be  derived  from  it.  For  example,  Judd  (1967) 
lists  nine  different  kinds  of  reflectance  functions  based  on  the  angular 
extent  of  the  incident  and  reflected  radiation.  Ail  of  them  can  be 
derived  from  the  BRDF. 

Note  that  the  BRDF  is  a four- dimensional  quantity  that  can  be 
thought  of  as  an  infinite  family  of  two-dimensional  light  distribution 
functions- -one  for  every  possible  angle  at  which  the  incident  beam  can 
strike  the  surface  element.  This  involves  an  overwhelming  quantity  of 
data,  especially  where  high  directional  resolution  is  needed  to 
describe  glints  and  specularities. 
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